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ABSTRACT 

We present new constraints on the relationship between galaxies and their host dark 
matter halos, measured from the location of the peak of the stellar-to-halo mass ra¬ 
tio (SHMR), up to the most massive galaxy clusters at redshift z ^ 0.8 and over 
a volume of nearly 0.1 Gpc^. We use a unique combination of deep observations in 
the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplemented by 
60 000 secure spectroscopic redshifts, analysing galaxy clustering, galaxy-galaxy 
lensing and the stellar mass function. We interpret our measurements within the halo 
occupation distribution (HOD) framework, separating the contributions from cen¬ 
tral and satellite galaxies. We find that the SHMR for the central galaxies peaks at 
-Wh.peak = l-9*o’i ^ IO^^Mq with an amplitude of 0.025, which decreases to ~ 0.001 
for massive halos (Mh > IO^^Mq). Compared to central galaxies only, the total SHMR 
(including satellites) is boosted by a factor 10 in the high-mass regime (cluster-size 
halos), a result consistent with cluster analyses from the literature based on fully in¬ 
dependent methods. After properly accounting for differences in modelling, we have 
compared our results with a large number of results from the literature up to z = 1: 
we find good general agreement, independently of the method used, within the typical 
stellar-mass systematic errors at low to intermediate mass (M* < IO^^Mq) and the 
statistical errors above. We have also compared our SHMR results to semi-analytic 
simulations and found that the SHMR is tilted compared to our measurements in such 
a way that they over- (under-) predict star formation efficiency in central (satellite) 
galaxies. 

Key words: cosmology: observations - dark matter - galaxies: clusters: general - 
gravitational lensing: weak. 


1 INTRODUCTION 

The last few years have seen an increasing interest in sta¬ 
tistical methods linking observed galaxy properties to their 
dark matter halos, owing to the availability of numerous 
large scale multi-wavelength surveys. Those techniques are 
based on the assumption that the spatial distribution of dark 
matter is predictable and one is able to match its statistical 
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properties with those of the galaxies. The halo model (see 
ICoorav fc Shethl2002ll is a quantitative representation of the 
distribution of dark matter, characterised by three main in¬ 
gredients: the halo mass function describing the number den¬ 
sity of halos per mass, the halo bias tracing the clustering 
amplitude, and the halo density profile. 

Galaxies are born and evolve in individual halos where 
the baryonic gas condensates, cools and forms stars. Galax¬ 
ies are gravitationally bound to dark matter and share a 
common fate with their host, e.g. during mergers. Although 
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we understand qualitatively individual physical processes 
likely to be involved in galaxy evolution, a number of key 
answers are missing. 

Observations show that a fraction of galaxies experi¬ 
enced star formation quenching and have become passive, 
shaping t he galaxy populat i on into a bimoda l blue/red dis¬ 
tribution dFaber et al.ll200^ : lllbert et al.ll2013l i. The number 
of these passive galaxies is higher today than in the past and 
increases with increasing halo mass. Might feedback pro¬ 
cesses in massive halos be responsible for this, or is there 
a universal critical stellar mass above which star formation 
ceases, independently of the halo mass? Studying the con¬ 
nection between galaxies and their host halos is crucial to 
answer these questions. 

Another enigmatic question is the low stellar mass 
fraction in low mass halos, se en in early studies connect¬ 
ing g a laxies to their host halos llYang, Mo fc van den Boschl 
12003 : IVale fc Ostrikerl l2006l : IZheng. Coil fc Zehavil l2007li . 
In fact, when measuring the stellar-to-halo mass ratio 
(SHMR) as a function of time, we observe that stellar 
mass is building up asymme trically, first in massive halos, 
later on in low-mass halo s dConrov. Wechsler fc KravtsovI 
l2006l : [Behroozi et al.ll20l3 l. This asymmetry in the SHMR 
is one corollary o f the so-called galaxy downsizing effect 
l|Cowie et al.lllOOfil i. In low-mass halos, stellar winds and su¬ 
pernovae may slow down star formation until the potential 
well grows deep enough to retain the gas and increase the 
star formation rate. Again, it becomes necessary to relate 
galaxy properties to their host halo mass. 

A number of studies have related galaxy properties to 
dark matter halos using the abundance matching technique 
ijMarinoni fc Hudsonl 20021: Conroy. Wechsler fc KravtsovI 
200^ Behroozi. Conroy fc Wechsler 2010l : Guo et ^ 2010l : 


Moster et al I 2 OIOI) , which employs the stellar mass (or lumi 


nosity) function and the halo mass function to match halo- 
galaxy properties based on their cumulative abundances. 
The conditional luminosity functi on technique proposed by 
lYang, Mo fc van den Boschl d2003li includes a parameterised 
M* — Mh relationship whose parameters are fitted to the 
luminosity function. Both this formalism and recent abun¬ 
dance matching studies feature a scatter in M* at hxed Mh, 
which is an important ingredient to account for, given the 
steep relation between the two quantities at high mass. 

More recently, models adopting a similar approach to 
abundance matching consist of directly populating dark 
matter halos from N-body simulations, to reproduce the 
observed stellar mass functions as a function of redshift, 
using a parameterised star formati on rate (SFR) model 
to account for redshift ev olution (iBehroozi et al.l l2013l : 
iMoster. Naab fc Whitel[2012li . 

Except in some rare cases where c entral or satellite 
galaxies can be ind ividually identihed fe.g. lMore et al.ll201ll : 
ICeorge et al.l[201lh . in studies based on luminosity or stel¬ 
lar mass distributions, the satellite galaxies’ properties can¬ 
not be disentangled from those of the central galaxies. To 
remedy the problem, abundance matching techniques either 
assume an ad-hoc fraction of satellites or use a sub-halo 
mass function estimated from numerical simulations. Unfor¬ 
tunately, as sub-halos may be stripped and disappear after 
being accreted onto larger halos, the subhalo mass func¬ 
tion at the time considered might not correspond to the 
distribution of satellites, and one must consider the mass of 


sub-halos at the time of accretion, further extrapolated to 
the time considered. Obviously these complications limit the 
amount of information one can extract about galaxy satel¬ 
lites. 

Galaxy clustering, on the other hand, allows separation 
of the contributions from central and satellite galaxies due to 
the different typical clustering scales. To model the cluster¬ 
ing signal of a given galaxy population, the halo occupation 
distribution (HOD) formalism assumes that the galaxy num¬ 
ber per halo is solely a function of halo mass and that the 
galaxy s atellite distribution is correlated to that of the d ark 
matter llBerlind fc Weinberdl20o3 : iKravtsov et ^120041 '). 

One achievem ent of HOD modelling was to demonstrate 
from simulations llBerlind et al.l l2003l : iMoster et al.l l201Clh 
that only a handful of parameters was necessary to fully 
describe galaxy halo occupation. This parametric HOD was 
fitted to a number of observations over a large range of red- 
shifts and galaxy properties. Among the more remarkable 
results are the local Universe galaxy clustering and abun¬ 
dance matching studies performed on the Sloan Digital Sky 
Survey (see e.g . Zehavi_et_^ 201 ill and at higher redshifts 
llFoucaud et ai.ll2010l: Wake et al. 201ll: ICoupon et al.ll2012l : 


Ide la Torre et al~ 2013l~ Martinez-Manso et al. 20l4l. 


However, some underlying assumptions on the distribu¬ 
tion of dark matter halos implied in the HOD formalism are 
observationally challenging to conhrm and one has to rely 
on N-body simulations. Fortunately, additional techniques 
may be used to relate galaxy properties to halo masses, 
among which gravitational lensing is one of the most pow¬ 
erful probes: by evaluating the distortion and magnification 
of background sources, one is able to perform a direct es- 
timation of the dark matter h alo profile (for a review, see 
iBartelmann fc Schneideill200ll l. The low signal-to-noise ra¬ 
tio associated with individual galaxies, however, forces us 
to “stack” them (e.g. binned together within narrow stellar 


ina jBrainerd. Blandford & Smai 

19961: Hudson et al.l 

19981; 

Hoekstra. Yee & GladdersI l2004l: 

Mandelbaum et al.l 

2005I; 

Yoo et al.ll2006l: van Uitert et al. 

201l|: 

Velander et al.l 

2014|; 

Cacciato. van Uitert & Hoekstra 

2014] 

Hudson et al. 20151'). 


Clearly, each of the above methods brings a differ¬ 
ent piece of information and combining all observables to¬ 
gether is particularly interesting, although doing so prop- 
erly is challengin g . In a recent study using COSMOS data, 
iLeauthaud et all ll2012h have successfully combined galaxy 
clustering, galaxy-galaxy lensing and the stellar mass func¬ 
tion, fitted jointly and interpreted within the HOD frame¬ 
work: the authors have used a global central galaxy M* — Mh 
relationship (as opposed to measuring the mean Mh per bin 
of stellar mass) and extended it in a consistent way to satel¬ 
lite galaxies. 

In this paper, we apply this advanced formalism using 
a new dataset covering a uniquely large area of ~25 deg^ 
with accurate photometric redshifts in the redshift range 
0.5 < 2 < 1 and stellar masses > lO^^M©. Our galaxy 
properties’ measurements are calibrated and tested with 
70 000 spectroscopic redshifts from the VIPERS survey and 
a number of publicly available datasets. Our data span a 
wide wavelength range of ultra-violet (UV) deep data from 
GALEX, optical data from the CEHT Legacy Survey, and 
As-band observations with the CEHT WIRCam instrument. 
This large statistical sample allows us to measure with high 
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precision the stellar mass function, the galaxy clustering, 
and we use the CFHTLenS shear catalogue to measure 
galaxy-galaxy lensing signals. The galaxy clustering is mea¬ 
sured on the projected sky for the photometric sample and 
in real space for the spectroscopic sample. 

This paper is organised as follows: in Section [51 we de¬ 
scribe the observations, the photometric redshift and stel¬ 
lar mass estimates. In Section [3] we present the measure¬ 
ments of the stellar mass function, the galaxy clustering 
(both from the photometric and spectroscopic samples) and 
galaxy-galaxy lensing signals. In Section [4] we describe the 
HOD model, and the MCMC model fitting results are given 
in Section (5] In Section [6l we discuss our results and con¬ 
clude. Throughout the paper we adopt the following cosmol- 
ogy: Ho—72 krn s“^ Mpc“^ and flm = 0.258, Da = 0.742 
(iHinshaw et al.l[2009l ) unless otherwise stated. To compute 
stellar m a sses w e adopt the initial mass function (IMF) of 
IChabrieil (120031 ') truncated at 0.1 and l OOMf?), and the stel¬ 
lar po pulation synthesis (SPS) models of iBruzual &: Chariot! 
ll2003ll . All magnitudes are given in the AB system. The dark 
matter halo masses are denoted as Mh and defined within 
the virial radius enclosing a mean overdensity Avir compared 
to the mean density backg r ound, taking the formula from 
IWeinberg fc Kamionkowskj (120021 '). At a = 0.8, Avir = 215. 
All masses are expressed in unit of M©. Measured quantities 
are denoted as X and theoretical quantities as X. We call 
cosmic variance the statistical uncertainties caused by the 
density fluctuations of dark matter and we define the sam¬ 
ple variance as the sum of the cosmic variance and Poisson 
noise variance. 


2 DATA 

In this work, we combine several datasets to build a volume- 
limited sample of galaxies more massive than M* = 10^° M© 
in the redshift range 0.5 < 2 < 1. Our galaxy selec¬ 
tion is based on NIR {Ks < 22) observations, collected 
in the two fields of the VIMOS Public Extragalactic Red¬ 
shift Survey (“VIPERS-Wl” and “VIPERS-W4”), overlap¬ 
ping the CFHTLS-Wide imaging survey, and covering a to¬ 
tal unma sked area of 23.1 deg^. We refer to lArnouts et al] 
(!in preD.h for a complete description of the multi-wavelength 
UV and NIR observations, reduction and photometry. 

Our background galaxy selection used for the measure¬ 
ment of the lensing signal is based on the CFHTLS-Wide 
i-band selection in the area that overlaps with the NIR ob¬ 
servations. 


2.1 The CFHTLS-Wide survey 

The Canada-France-Hawaii Telescope Legacy Surve3{3 
(CFH TLS) is a photometr ic survey performed with Mega- 
Cam (iBoulade et al.l[2003l l on the CFHT telescope in five 
optical bands u*, g,r,i, z {i < 24.5 — 25, 5cr detection 
in 2" apertures) and covering four independent patches 
in the sky over a total area of 154 deg^. In this anal¬ 
ysis, we use the photometric and shear catalogues pro - 
duced by the CFHTLen^ team (iHevmans et al.l l2012ll . 


The CFHTLenS photom etry is performed with SExtractor 
( Bertin fc Arnorit3 19961') on the PSF-hom ogenised images 
( Erben et al.l2013l : lHildebrandt et al.l2012l ). Magnitudes are 
based on the MAG_IS0 estimator where the isophotal aper¬ 
tures are derived from the i-band detection image. This 
approach optimises the colour measurements and leads 
to an improvement in t he photometric redshift accuracy 
dHildebrandt et irDl2012ll . To estimate the total magnitude 
of each source, a global shift is applied to the MAG_ISD mag¬ 
nitude in all the bands based on the difference between 
MAG_IS0 and MAG_AUTQ magnitudes, as measur ed in the i- 
band detection image (iHildebrandt et ^l2012ll . 

As the magnitude errors are measured with SExtractor 
directly from the local background in the PSF-homogenised 
image, we need to correct for the noise correlation intro¬ 
duced by the convolution process. To do so we multiply 
the CEHTLenS magnitude errors in all bands by the ra¬ 
tio of the i-band detection image errors to the i-band PSF- 
homogenised image errors. The correction factor ranges from 
3 to 5, where the strongest correction occurs when the seeing 
difference between the i-band and the worse-seeing image is 
the largest. As the i-band image is usually the best-seeing 
image, this procedure may slightly overestimate the correc¬ 
tion in the other bands, however we neglect it here. 

In addition, magnitude errors must be rescaled to ac¬ 
count for image resampling. Two independent tests have 
been performed to accurately estimate the correction fac¬ 
tor: we measured the dispersion of magnitudes between 
the i-band detection (un-convolved) magnitudes and the 
CFHTLS-Deep magnitudes, and between duplicated obser¬ 
vations of the same object in the overlapping regions of ad¬ 
jacent tiles. We find that the errors must be rescaled by a 
factor of 2.5. 

The footprints of the CFHTLS MegaCam tiles overlap¬ 
ping the VIPERS survey are shown as gray squares in Eig. [T] 


2.2 The Near-IR observations 

We have conducted a As-band follow-up of t he VIPERS 
fields with the WIRCam instrument at CFHT (jPueet et al.l 
l2004f) for a total allocation time of ~120 hours. The integra¬ 
tion time per pixel was 1050 seconds and the average see¬ 
ing of all the individual exposures was 0.6"±0.09. The data 
have been reduced by the Terapix tearr0: the images were 
stacked and res ampled on the pixel grid of the CFHTLS- 


T0007 release ([Hudelot et al. 

2012 !). The images reach a 

depth of As = 22 at ~ 3(7 ( 

Arnouts et aLllin nreD.!). The 


photometry was performed wit h SExtractor in dual image 
mode with a gri — image dSzalav. Connolly fc Szokolvl 
I 1999 II as the detection image. To correct for the noise cor¬ 
relation introduced by image resampling, we multiply the 
errors by a factor 1.5, obtained from the dispersion be¬ 
tween the WIRCam As-band magnitudes and the magni¬ 
tudes measured in t he deeper (A < 24.5) UKIDSS Ultra 
Deep Survey fUDS. [Lawrence et al.l [2()()n L We also used 
the UDS survey to confirm that our sample completeness 
based on gri — detections reaches 80% at Ks = 22. Us¬ 
ing the WIRCAM/CFHTLS-Deep data with an f-band cut 
simulating the CFHTLS-Wide data depth, we have checked 


^ http://www.cfht.hawaii.edu/Science/CFHTLS/ 

^ http://cfhtlens.org/ ® http://terapix.iap.fr/ 
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Figure 1. Footprints of the different datasets used in this work. Our selection is based on WIRCam data shown in red and covering 
approximately 25 deg^ (23.1 deg^ after masking). The CFHTLS MegaCam pointings are shown in grey, the GALEX DIS observations 
as large blue circles (in purple if overlapped with WIRCam), the spectroscopic surveys VIPERS/VVDS in light green and PRIMUS in 
dark green. The SDSS/BOSS coverage is almost complete. The data outside the WIRCam footprint are not used, and shown here only 
for reference. 
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that this incompleteness is caused by red galaxies above 
z = 1 and does not affect our sample selected in the range 
0.5 < z < 1. The Ks MAG_AUTD estimates are then simply 
matched to their optical counterparts based on position. 

In addition t o this dataset, we also use the CFHTLS- 
D1 WIRDS data l|Bielbv et al.lEoi^ . a deep patch of 0.49 
deg^ observed with WIRCam J-, H- and As-bands and cen¬ 
tered on 02^26™59'’, —04°30'00”. All three bands reach 50% 
completeness at AB magnitude 24.5. 

The WIRCam observations are shown in Fig. [T] as the 
red regions. After rejecting areas with poor WIRCam pho¬ 
tometry and those with CFHTLenS mask flag larger than 2, 
the corresponding effective area used in this work spans over 
23.1 deg^, divided into 15 and 8.1 deg^ in the VIPERS-Wl 
and VIPERS-W4 fields, respectively. 

2.3 The UV-GALEX observations 

When available, we make use of the UV deep imaging 
photometry from the GALEX satellite llMartin et al.l [2005c 
iMorrissev et al.ll2005l ). We only consider the observations 
from the Deep Imaging Survey (DIS), which are shown in 
Fig. [T]as blue circles (0 ~ 1.1°). All the GALEX pointings 


were observed with the NUV channel with exposure times 
of Texp ^ 30 ksec. FUV observations are available for 10 
pointings in the central part of Wl. 

Due to the large PSF (FWHM~5"), source confusion 
becomes a major issue in the deep survey. To extract the UV 


photometry we use 

a dedicated photometric code, EMphot 

llConseil et al.l 2011 

) which will be described in a separate 

oaner (Vibert et al. 

lin DreD.h. In brief. EMnhot uses I/-band 

(here the CFHTLS w-band) detected objects as a prior on 


position and flux. The uncertainties on the flux account for 
the residual in the [simulated—observed] image. The images 
reach a depth of miyuv ~ 24.5 at ~ 5a. As for the WIRCAM 
data, the GALEX sources are matched to the optical coun¬ 
terparts based on position. 

The NUV observations cover only part of the WIRCam 
area with ~10.8 and 1.9 deg^ in VIPERS-Wl and VIPERS- 
W4, respectively. The UV photometry slightly improves the 
precision of photometric redshifts and the stellar mass es¬ 
timates in the GALEX area. However, by comparing our 
measurements inside and outside the GALEX area, we have 
checked that the addition of UV photometry does not make 
a significant change for the galaxies of interest in this study. 
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Therefore, in the final sample, we mix galaxies inside the 
GALEX area with those outside. 


2.4 Spectroscopic data 

To optimise the calibration and the validation of our pho¬ 
tometric redshifts we make use of all the spectroscopic red- 
shifts available in the WIRCam area. 

The large st sample is based on the VIPERS spectro- 
scopic survey iIGuzzo et al.ll2oi3 : iGarilli et al.ll20f3l and its 
first public data release PDRllj. VIPERS aims to measure 
redshift space distortions and explore massive galaxy prop¬ 
erties in the range 0.5 < a < 1.2. The survey is located in 
the W1 and W4 fields of the GFHTLS-Wide survey and 
will cover a total area of 24 deg^ when completed, with 
a sampling rate of ~40% down to i < 22.5. In Fig. [T] we 
show the layout of the VIMOS pointings as the light-green 
squares. The PDRl release includes redshifts for ~ 54 204 
objects. After keeping galaxy spectra within the WIRGam 
area (44 474) and with the highe st confidence flags between 
2.0 and 9.5 (95% confidence, see lGuzzo et al.ll2014l i. we are 
left with 35 211 galaxies, which corresponds to a spectro¬ 
scopic success rate of 80%. 

In addition to VIPERS, we also consider the following 
spectroscopic surveys: 

• the VIM OS-VLT Deep Survey (WP S) F02 and Ultra- 
Deep Survey dLe Fevre et al.lboOSl . lioi^ i which consist of 
11353 galaxies down to i < 24 (Deep) and 1125 galaxies 
down to i < 24.5 (Ultra-Deep) over a total area of 0.75 deg^ 
in the VIPERS-Wl field. We also use part of the VIMOS- 
VLT F22 Wide Su rvey with 12 995 galaxies over 4 deg^ 
down to i < 22.5 jGarilli et al.l l2008l . shown as the large 
green square in the southern part of the VIPERS-W4 field 
in Fig. [T]). In total, we use 5 122 galaxies with secure flags 3 
or 4 from the VVDS surve ys within the WI RCam area; 

• the PRIMUS survey JCoil et al.ll201l| i which consists 
of low resolution spectra (A/A A ~ 40) for galaxies down 
to i ~ 23 and overlapping our VIPERS-Wl field. PRIMUS 
pointings are shown as the dark green circles in Fig. [T] We 
keep 21 365 galaxies with secure flags 3 or 4; 

• the SD SS-BOSS spectro scopic survey based on data re¬ 
lease DRIO (lAhn et alJl2014h down to i < 19.9, overlapping 
both VIPERS-Wl and VIPERS-W4 helds, totalling 4 675 
galaxies with zWarning=0 (99% conhdence redshift) within 
our WIRCam area. 

In total, the spectroscopic sample built for this study 
comprises 62 220 unique galaxy spectroscopic redshifts with 
the highest confidence flag. We use the spectroscopic red¬ 
shift value, when available, instead of the photometric red¬ 
shift value. The galaxies with a spectroscopic redshift repre¬ 
sent 6.5% of the full sample, and 12% after selection in the 
range 0.5 < z < 1, where most of the galaxies are from the 
VIPERS sample. 


Table 1. Magnitude zero-point offsets measured per CFHTLS 
MegaCam pointing in VIPERS-Wl and VIPERS-W4 (mean and 
standard deviation). J and i7-band zero-points were computed 
for the pointings overlapping WIRDS data. 


Filter 

VIPERS-Wl 

VIPERS-W4 

FUV 

0.18 ± 0.11 

0.02 ± 0.16 

NUV 

0.11 ± 0.09 

0.15 ± 0.10 

U 

0.10 ± 0.03 

0.13 ± 0.03 

g 

-0.02 ± 0.01 

-0.01 ± 0.01 

r 

0.02 ± 0.01 

0.01 ± 0.01 

i 

-0.01 ± 0.01 

-0.00 ± 0.01 

z 

-0.02 ± 0.01 

-0.01 ± 0.01 

j 

0.08 ± 0.05 

— 

H 

0.02 ± 0.05 

— 

K 

0.02 ± 0.03 

0.01 ± 0.05 


2.5 Photometric redshifts 


To compute the photome tric redshifts, we use the tem¬ 
plate fitting code LePhar43 (lArnouts et al.l 19991 : lllbert et al.l 


20061). We adopt s imilar extinction laws and parameters as 


jlbert et al.l (l2009l) used in the C OSMOS field ijs coville et al.l 
2007h . and identical priors as in ICoupon et al.r i 2009h based 


on the VVDS redshift distribution and maximum allowed g- 
band absolute magnitude. We note that the use of priors is 
essential for the 2 > 1, low signal-to-noise (or no NIR flux), 
galaxies used as lensed (background) sources (see Section|3). 
A probability distribution function (PDF) in steps of 0.04 
in redshift is computed for every galaxy. 

We use the full spectroscopic sample to adjust the mag¬ 
nitude relative zero-points in all the passbands on a Mega- 
Cam pointing-to-pointing basis. For the pointings with no 
spectroscopic information, we apply a mean correction ob¬ 
tained from all the pointings with spectra. The mean zero- 
point offsets and standard deviations in all passbands are 
given in Table [T] for the two helds separately. We further add 
the zero-point scatter in quadrature to the magnitude errors 
in each band. We recall that these zero-point corrections may 
not represent absolute calibration offsets but rather relative 
(i.e. depending on colours) ones and tied to the adopted 
spectral energy distribution (SED) template set. We come 
back to the impact of this issue on stellar mass measure¬ 
ments in Section m 

Our SED tem plates are based on the library used in 
lllbert et al.l ll2009l) . however the fewer bands used in this 
study compared to COSMOS necessitate adapting the tem¬ 
plates to reduce redshift-dependent biases. The init i al tem - 
plates are based on the SEDs from IPolletta et al.l ll2007l) , 
complemented by a num ber of starburst SEDs from the 
iBruzual &: Charlotl (l2003l) SPS library. Using 35 211 spec¬ 
troscopic redshifts from VIPERS, we adapt the templates 
with LePhare using the following procedure. First, a best-ht 
template from the original set is found for each galaxy and 
normalised to unity, and the photometry is then corrected 
into the rest frame given the spectroscopic redshift value. 
The rest-frame photometry for all galaxies with identical 
best-fit templates is combined and the adapted template 
is constructed from the sliding-window median values as a 


http://vipers.inaf.it/rel-pdrl.html 


® http://www.cfht.hawaii.edu/ arnouts/lephare.html 
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function of wavelength. The process is repeated iteratively. 
Given the high number of galaxies with spectroscopic red- 
shifts, we found that only two iterations were necessary to 
reach convergence. Interestingly, although the improvement 
in the photometric redshift bias is significant, the new tem¬ 
plates appear very similar “by eye” compared to the original 
ones, which implies that small features in the SED tem- 
pl ates may lea d to la rge photometric errors, as also noted 
bv lllbert et al.l l|200fih . 

In Fig. (2] we show the accuracy of the photometric red- 
shifts by comparing with the spectroscopic redshift sample 
from VIPERS (i < 22.5, left panel) and VVDS Deep/Ultra- 
Deep (22.5 < i < 24.5, right panel). We observe a disper- 
sioqjof (j/(l-|-z) ~ 0.03 — 0.04 andafraction of catastrophic 
redshifts (| A^l ^ 0.15(1 -I- z)) of ?7 ~ 1 — 4%. The dispersion 
in both magnitude ranges is signif icantly better than p revi- 
ous results in the CFHTLS-Wide dCoupon et al.l[2009ll . due 
to th e choice of isophotal mag nitudes and PSF homogenisa¬ 
tion llHildebrandt et al.ll2012ll at faint magnitude, and the 
contribution of NIR data above 2 ~ 1. We note that the faint 
sample is compared to the VVDS redshifts where deep NIR 
data from WIRDS are available over a small part (< 1 deg^) 
of the field, and with a magnitude distribution biased to¬ 
wards bright galaxies compared to the photometric sample. 
Therefore, we foresee degraded photometric redshift perfor¬ 
mance elsewhere, mainly relevant for z > 1 galaxies. How¬ 
ever, as shown in Appendix El no systematic bias affecting 
our lensing measurements is introduced by the use of sources 
beyond z = 1. 


2.6 Stellar mass estimates 

To compute stell a r mas ses we adopt the same procedure 
as lArnouts et al.l (l2013l l and described in detail in their 
Appendix A. In brief, we use the photometric or spectro¬ 
scopic (when available) redshift and perform a mini- 
mization on a SED libra ry based on the SPS code from 
iBruzual fc Charlotl (l200.3l l. The star formation history is ei¬ 
ther constant or described with an exponentially declining 
function, with e-folding time 0.01 Si r ^ 15. We use two 
metallicities {Zq, 0.2 Z(7)) and adop t the IChabrierl ll2003l f 
IMF. As discussed in lArnouts et al.l (|2013f l. the use of var¬ 
ious dust extinction laws is critical to derive robust SFR 
and stellar mass, and in this work we adop t their choices 
for d iffering attenua tion curves: a starb urst dCalzetti et al.l 
I 2 OO 0 I I. a SMC-like (IPrevot et al.l [iQS^l and an intermedi¬ 
ate slope (A~° ®) law. We consider reddening excess in the 
range 0 ^ E{B — V) ^ 0.5. When fixing the redshift, the 
typical 68% stellar mass statistical uncertainty, as derived 
by marginalising the likelihood distribution, ranges from 
cr(M*) ~ 0.05 to 0.15 for galaxies with Ks ^ 22 and z < 1. 
This stellar mass uncertainty is an underestimate, since we 
neglect photometric redshift uncertaintiefl. 

In addition to statistical errors, in Section 13.51 we in¬ 
vestigate the different sources of systematic effects in the 


^ Defined as the normalised median absolute deviation 
dHoaglin. Mosteller fc TukevI [l983h : 1.48XMedian)l^s — Zp\/(\- + 
^s)), and robust to outliers. 

^ We will see in Section l4.2l that our model accounts for such an 
extra source of uncertainty in stellar mass through a stellar-mass 
dependent parameterisation of the stellar-to-halo mass scatter. 


stellar mass estimates, arising from our lack of knowledge 
of galaxy formation and evolution. The choice of differing 
dust trea tments (and r esulti ng dust attenuation laws) is one 
of them: lllbert et al.l (l2010l l have measured a shift of 0.14 
dex, with a large scatter , betwe en stellar masses estimated 
with th e ICharl ot &: Falll ll2000l l dust prescription and the 
ICalzetti et al. (2000 ll aHenuation law. The dust parameter¬ 
isation leads to systematics larger than the statistical errors 
in the stellar mass function. Even more critical is the choice 
of the SPS model and the IMF (see more detailed system¬ 
atic errors analysis in fBehroozi, Conroy fc Wechsle: 


s ystem- 

dH; 


iMarchesini et alll2009l : iFritz et al.ll2014jl ' leading to svstem- 
atic differences in stellar mass estimates up to 0.2 dex. One 
must keep these limitations in mind when comparing results 
from various authors using different methods, and we come 
back to these issues when presenting our results. 


3 MEASUREMENTS 

We aim to compute high signal-to-noise measurements of 
four distinct observables: the stellar mass function <))(M*), 
the projected galaxy clustering w{9), the real-space galaxy 
clustering Wp(rp), and the galaxy-galaxy lensing AE(r). 

To do so, we select volume-limited samples in the red¬ 
shift range 0.5 < z < 1, where the high sampling rate of 
VIPERS and our NIR data guarantee both robust photo¬ 
metric redshift and stellar mass estimates. As for the stellar 
mass function, we adopt a lower mass limit of M* = 10^° M© 
and employ the Umax estimator to correct for galaxy incom¬ 
pleteness near z = 1. The total volume probed in this study 
is 0.06 Gpc®. 

The stellar mass bins for the clustering and lensing 
measurements are defined to keep approximately a constant 
signal-to-noise ratio across the full mass range (which may 
lead to differing mass cuts depending on the observable), 
and guarantee complete galaxy samples (see Appendix lAl). 
We summarise our samples’ properties in Tabled] 

To measure each of the observables described below, we 
use the parallelised code SWOT, a fast tree-code for computing 
two-point correlations, histogra ms, and galaxy-galaxy lens¬ 
ing signals from large datasets (ICoupon et al.ll2012I) . The 
stellar mass function is expressed in comoving units, whereas 
the clustering and galaxy-galaxy lensing signal are measured 
in physical units. We estimate statistical covariance matri¬ 
ces from a jackknife resampling of 64 sub-regions with equal 
area (0.35 deg^ each), by omitting a sub-sample at a time 
and computing the properly normalised standard deviation 
(see more details in ICoupon et al.ll2012h . This number was 
chosen to meet both requirements of using large enough sub- 
regions to capture the statistical variations at large scale, 
while keeping a sufficient number of sub-samples to compute 
a robust covariance matrix. Nevertheless, we expect the pro¬ 
jected galaxy clustering errors to be slightly underestimated 
on scales larger than the size of our sub-regions, ~ 0.5 deg, 
and the noise in the covariance matrix to potentially bias 
the best-fit value. 

A random sample with 1 million objects is constructed 
using our WIRCam observations layout and the union of 
the WIRCam and CFHTLenS photometric masks. For real- 
space clustering, measured from VIPERS spectroscopic red¬ 
shifts, the random sample is constructed using the layout of 
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Figure 2. Photometric redshifts measured with ugrizK (left) or ugrizJHK (right) photometry versus VIPERS and VVDS spectroscopic 
redshifts. Left: 17.5 < i < 22.5, where the sample is dominated by galaxies between 0.5 < 2 : < 1.2 due to the VIPERS selection. Right: 
22.5 < i < 24.5, from the VVDS Deep and Ultra-Deep surveys. The limits for the outliers are shown as red dotted lines. 


Table 2. Sample mass definitions in log(M*/M 0 ) and number of galaxies in each sample. The parent sample comprises a total of 352 585 
galaxies. 


Sample 

Clustering - 
Mass cut 

- 'w{6) 
Number 

Clustering - 
Mass cut 

Wp{rp) 

Number 

Lensing 
Mass cut 

Number 

1 

10.00 - 10.40^ 

23 886 

10.60 - 10.90^ 

2154 

10.00 

- 10.401 

23 886 

2 

10.40 - 10.60 

36 560 

10.90 - 11.20^ 

1964 

10.40 

- 10.65 

45 032 

3 

10.60 - 10.80 

31900 

11.20 - 12.00 

816 

10.65 

- 10.80 

23 427 

4 

10.80 - 11.00 

24451 

— 

— 

10.80 

- 10.95 

19 293 

5 

11.00 - 11.20 

13 538 

— 

— 

10.95 

- 11.15 

16317 

6 

11.20 - 12.00 

6 326 

— 

— 

11.15 

- 12.00 

8 654 


10.5 < 2 < 0.7 
2 0.5 < 2 < 0.8 


the VIPERS PDRl geometry (and photometric masks) plus 
a random redshift drawn in the range 0.5 < 2 < 1 from a 
distribution following dV/dz, to match our volume limited 
samples. The sub-regions for the measurements of statistical 
errors are constructed by SWOT based on the random cata¬ 
logue: the field is divided into 64 areas with an equal number 
of random objects. 


3.1 The stellar mass function 

The stellar mass function (j){M0 = dn/dlog M^, is mea¬ 
sured per unit of comoving volume in 10 equally spaced 
logarithmic mass bins of width 0.2 dex, centered on the 
mass mean weighted by the number of galaxies. To cor¬ 
rect for the incompleteness in the low-mass galaxy sample 
(10^° < M*/M0 < 10'^° "') occurring near 2 = 1 (see Ap¬ 
pendix we up-weight low-redshift galaxies by a factor 
1/Vjnax defined as: 


Rm 



—d2, 

d2: 


( 1 ) 


where fl is the solid angle of the survey, 23.1deg^, V the co¬ 
moving volume per unit area, and 2max the maximum red- 
shift for a galaxy to be observed given a Ks < 22 magnitude 
cut, calculated with LePhare. 

We have performed a number of tests to check our in¬ 
ternal error estimates. In the top panel of Fig.|3]we show our 
stellar mass function error estimates (square root of the co- 
variance matrix diagonal) as a funct ion of stellar mass co m- 
pared to the getcv code estimate of iMoster et al.l (1201 il l at 
2 = 0.8. The latter code computes the theoretical expec¬ 
tations of cosmic varianc^l assuming a prediction for dar k 
matter clustering and galaxy biasing (iBardeen et al.lHo^ l. 
We add to the getcv cosmic variance the theoretical Pois¬ 
son error and show the resulting (total) sample variance as 
the thick line in the bottom panel. Our jackknife estimate 


® We note that the highest mass bin galaxy bias was estimated a 
posteriori from our HOD results, since it was not provided by the 
authors of getcv, although the contribution of cosmic variance is 
negligible compared to the Poisson error in this bin, populated 
by rare massive galaxies. 
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Figure 3. Stellar mass function statistical errors as function of 
stellar mass (top) and area (bottom). In the top panel we show 
the jackknife estimator based on 64 sub-regions and multiplied 
by a factor 2, compared to the theoretical cosmic v ariance plus 
Poisson error derived from the iMoster et alj (l201lh getcv code 
(the Poisson error only is shown as the dotted line). The bottom 
panel shows an alternative internal estimate based on the stan¬ 
dard deviation of sub-regions as a function of their size, in two 
mass bins (logM*/MQ = 10.10 and 11.89), extrapolated to the 
size of the full survey (dashed lines in both panels). As in the top 
panel, the black dots are the jackknife estimates, for which the 
cosmic variance part has been multiplied by 2. 


is represented as the black points, for which we find that 
the cosmic variance part (after subtracting Poisson noise) 
needs to be multiplied by a factor 2 to agree with theoreti¬ 
cal expectations (we then multiply the covariance matrix by 
a factor 4). We have not found a definitive explanation for 
the underestimation of the errors from the jackknife resam¬ 
pling, however it is likely caused by the strong correlation 
between bins (a combined effect of stellar mass scatter and 
large scale structure correlations). 

In the bottom panel of Fig. [S] we show an alternative 
internal estimator as function of area, based on the standard 
deviation of subsamples with sizes varying from 0.1 deg^ to 
2.9 deg^ (the black dots represent our Jackknife estimates in 
the two mass bins {log /Mq) = 10.10 and 11.89). We use 
a power law fit (the amplitudes of the error-bars are arbitrar¬ 
ily scaled to the square root of the number of subsamples, 
ranging from %/256 to %/&) to extrapolate to the full size of 
the survey. The extrapolated values are shown as the dashed 
line in the top panel of Fig. [S] The bin correlations between 
small sub-samples may tilt the slope of the fit and lead to an 
overestimate of the extrapolated error estimate, as observed 
in the low mass bin. In the high mass bins, characterised 


by an uncorrelated sampling variance dominated by Poisson 
noise, the extrapolated estimate is consistent with both the 
jackknife estimate and the theoretical Poisson noise. 


3.2 Projected galaxy clustering 

We measure the two-point correlation function w{6) in 10 
logarithmically spaced bins centered on the pair-number 
weighted averaged separation over the range 0.002° < 9 < 
2°. The modelled w{9) is compared to the measured 'w{9) 
by projecting the theoretical spatial clustering ^(r) onto the 
sample redshift distribution computed as the sum of photo¬ 
metric redshift P DFs (see Section 12.511. 

We use the iLandv fc Szalavl (1 19931) estimator follow- 
in g a similar pr o cedur e to that described in Section 3.3 
of ICoupon et ahl ll2012l) . Owing to the limited size of the 
survey, our measurements are affected by the integral con¬ 
straint, an effect that biases the clustering signal low. Here, 
we adopt a refined way to correct for it: the correction is cal¬ 
culated directly for every parameter set from the modelled 
'w{9) (instead of a pre-determined power law) a nd integrated 
over t he survey area using random pairs as in iRoche et al.1 
ll2002l) . leading to better agreement between the data and 
the model at large scales. Here the typical values of the in¬ 
tegral constraint range from 10“^ to 3 x 10“®. 

We have checke d, using the galaxy moc ks prepared for 
the VIPERS sample dde la Torre et al]l2013l) . that our jack¬ 
knife error estimates could reproduce within 20% the cor¬ 
rect sample variance amplitude of w{9) (this result is in 
agreement with a number of tests from the literature, e.g. 
IZehavi et al.l[2005l : iNorberg et al.l[2009l ). and we do not ap¬ 
ply any correction. 


3.3 Real-space galaxy clustering 

We measure the real-space galaxy clustering for the VIPERS 
spectroscopic sample by integrating the weighted redshift- 
space correlation function along the line of sight to alleviate 
redshift-space distortion effects: 

C(^P,phys,'7rphys)d7rphys , (2) 

where rp,phys and TTphys are the coordinates perpendicular 
and parallel to the line of sight, respectively, rp^phys is ex¬ 
pressed in physical coordinates and divided into 10 logarith¬ 
mically spaced bins centered on the pair-weighted averaged 
separation over the range 0.2 < rp,phys/Mpc < 10, and TTphys 
is divided into linear bins up to Timax = 40 Mpc. The value 
of TTmax is consisteutlj^ used in the derivation of the mod- 
elled Wp. As for w (9), ^(rp,phys, TTphys) is computed using the 
iLandv fc Szalavl estimator and the covariance matrix esti¬ 
mated from the jackknife resampling of 64 sub-regions. 

Each galaxy is weighted to account for the undersam¬ 
pling of the spectroscopic sample: we use the global colour 
sampling rate (CSR), target sampling rate (TSR), and suc¬ 
cess s ampling rate (SSR), as described in IPavidzon et ^ 
ll2013l) , to account for the VIPERS colour selection, the 
sparse target selection and measurement success as function 
of signal-to-noise ratio, respectively. In addition, we also use 
number-count normalised (to prevent global CSR, TSR and 
SSR double weighting) spatial weights computed for each 
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VIPERS panel bv ide la Torre et al.l ll2013h to correct for the 
position-dependent sampling. Here, the SSR is the most af¬ 
fected quantity, as a function of position in the sky, due to 
the differing observing conditions at the times of observa¬ 
tion. 

Small pair incompleteness due to “slit collision” is cor- 


rected by a factor 1 -I- wa, such that: 


1 , ~ _ 1 + rup 

-L “r 'i^p,corr ^ ? 

l + WA 

(3) 

where 


, ~ , 0.03 

1 -1- WA = 1- 

?"p,phys 

(4) 


is derived fro m the projected co r relati on as function of an¬ 
gular scale by Ide la Torre et all (l2013l i and translated into 
physical scales at 2 = 0.8. We note that given our conserva¬ 
tive small scale cut of rp^phys > 0.2, the correction remains 
below 15%. 


3.4 Galaxy-galaxy lensing 

The gravitational lensing signal produced by the foreground 
matter overdensity is quantified by the tangential distortion 
of background sources behind a sample of stacked “lens” 
galaxies, also k nown as the wei g hted galaxy-galaxy lensin g 
estimator (e.g. IYoo et akl I2OO6I : iMandelbaurn et S I2OO6II . 
The excess surface density of the projected dark matter halo 
relates to the measured tangential shear through: 

^E(rp^phys) — Ecrit X qt (rp^phys), (5) 

(see also Appendix We measure the signal in 10 log¬ 
arithmically spaced bins centered on the number-weighted 
averaged separation, in the range 0.02 < ^,phys/Mpc < 1. 
rp.phys is expressed in physical coordinateqj. 

The critical surface density Scrit is given by: 

V - -Pqs 

47rGN HoL Dls ’ ^ ’ 

with Dos the observer-source angular diameter distance, 
Dol the observer-lens (foreground galaxy) distance and Dls 
the lens-source distance. Gn is the gravitational constant 
and c the speed of light. All distances are computed in phys¬ 
ical coordinates using the photometric (spectroscopic when 
available) redshift. For photometric redshift values, a cut 
^source — ^lens > 0.1 X (1 -|- 2iens) is adopted. The background 
source galaxy sample includes all gala xies detected in th e 
i-band with a non-zero lensing weight (iMiller et al.ll20f^ . 
Here we do not restrict our redshift sample to Zp < 1.2, but 
consider galaxies at all redshifts, taking advantage of the 
improved photometric redshift estimates in our sample, in¬ 
creasing the background source sample by 30% compared to 
other CFHTLenS lensing studies, without introducing any 
systematic bias (see Appendix ICl). 

The galaxy shape measurement was performed on in- 
dividual exposure s using the lensfit analysis pipeline 
llMiller et al.l l2007l : iKitching et al.l l2008l : iMiller et al.l l2013l i 

® Note that the galaxy-galaxy lensing signal is measured in phys¬ 
ical units, whereas a number of authors assume comoving units, 
which would require multiplying the excess surface density by a 
factor of (1 + z)~‘^ compared to our definition. 


and sy stematics checks were conducted by iHevmans et al.l 
(l2012ll for cosmic shear (the projected large scale struc¬ 
ture lensing power spectrum). The lensing (inverse-variance) 
weights account for shape measurement uncert ainties 
dMiller et al.l l2013l ) . Following IVelander et al.l (120141 ') , who 
performed extensive systematics checks of the CFHTLenS 
shear catalogue specifically for galaxy-galaxy lensing (see 
their Appendix C), we do not reject those CFHTLS-Wide 
pointings that did not pass the requirements for cosmic 
shear, and we applied appropriate shape measurement cor¬ 
rections as described in their Section 3.1. 

We compute the boost factor (to account for dilu- 
tion due to sources physically associated wit h the lens, see 
ISheldon et al.l[2^004l : fMandelbaum et al.ll2006h by randomis¬ 
ing the source positions, and correct the final signal for it. 
On small scales (< 0.1 Mpc) the boost factor reaches up to 
20 % for the most massive galaxies. 

Here, the relatively low source density implies that our 
errors are dominated by the source galaxy shape noise, orig¬ 
inating from ellipticity measurement uncertainties and in¬ 
trinsic shape dispersion, rather than sample variance. In¬ 
deed, when compared to the sum of inverse-variance lensing 
weights, we have checked that our jackknife estimate was 
similar at all scales (with small off-diagonal correlation), 
confirming the negligible impact of cosmic variance (see Ap¬ 
pendix 0. 

Nevertheless, a correlation exists between the mass bins 
due to the re-use of background source galaxies. We neglect 
this contribution in the computation of the combined 
but we note that this correlation is likely to lead to under¬ 
estimation of our parameter confidence limits. 


3.5 Systematic errors in stellar mass 
measurements 


In this section we are concerned with systematic errors af¬ 
fecting the stellar mass measurements caused by the uncer¬ 
tainties in the assumed cosmology (i.e. volume and distance 
estimates), the dust modelling, and potential biases in the 
photometry. 

To assess the impact of systematics on the measure¬ 
ments of the observables, we propagate the errors affecting 
the stellar masses by changing one parameter configuration 
at a time, then re-computing all stellar masses and the ob¬ 
servables, and finally measuring the difference with the ref¬ 
erence measurements. We repeat the process for the three 
different kinds of systematics listed below: 


• assumed cosmology. We explore three ACDM parame¬ 
ter sets: in addition to the WMAP cosmology used in this 
study with Hp = 72 km s“^ Mpc“^, = 0.258, Ha = 0.742 
dHinshaw et al.l l2009h , a “concordance” cosmology model 
with Ho = 70 km s“^ Mpc“^, Hm = 0.3, Ha = 0.7, 
and the Planck cosmolog y with Hp = 67 km s~^ Mpc~^ , 
Hm = 0.320, Ha = 0.680 ( Planck Collaboration et al.ll20Lt ) 
are tested. In each case, the stellar masses and the observ¬ 
ables are consistently re-computed with the same cosmol¬ 
ogy. We note that the term “systematics” here refers to the 
choice for one or another set of parameters that produces a 
systematic shift in stellar mass and not to systematic errors 
associated to the measurement of cosmological parameters; 

• lens galaxy dust extinction modelling. We compute five 
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different stellar masses for each galaxy by varying one aspect 
at a time: two different extinction law configurations (among 
our choice of three laws, see Section 12.61) and three different 
E{B — V) maximum allowed values (ranging from 0.2 to 
0.7); 

• photometric calibration. As zero-point offsets do not 
correct for absolute calibration uncertainties (but do for 
colours), nor correct for photometric measurement biases 
(e.g. missing flux of bright objects), a change in the pho¬ 
tometric calibration may cause a shift in the best-fit tem¬ 
plate and further bias the stellar mass measurements. We 
re-compute stellar masses applying ad-hoc global shifts (in 
all bands) of —0.05 and -1-0.05 magnitude, which correspond 
to typical offsets caused by various calibration strat egies or 
photometry measurements l|Moutard et ahlfin DreD.l l. 

Results are shown in Fig. [H For each observable (top 
left: stellar mass function, top right: projected clustering, 
bottom left: real-space clustering, and bottom right: galaxy- 
galaxy leasing), we display the re-computed measurements 
divided by the reference quantities, in each of the “Cosmol¬ 
ogy”, “Extinction” and “Calibration” panels as well as the 
sum in quadrature of all these effects (“total”). The gray 
area corresponds to the maximum value among the differ¬ 
ing re-computations, not the standard deviation, as each of 
the solutions is equally likely to be opted for. Except for 
the stellar mass function, we only display the results in the 
most massive bins (where we observe the most significant 
changes), although the calculations were repeated in all mass 
bins. 

To allow comparison with the statistical errors, we over¬ 
plot the error bars from our jackknife error estimates. For 
the stellar mass function (whose jackknife error estimate is 
multiplied by a factor 2, see Section EIJ, the systematic er¬ 
rors compared to the statistical errors are striking, with the 
former being larger by one order of magnitude compared 
to the latter. The increase of the systematic errors towards 
the high mass regime is a direct consequence of the shift in 
stellar mass and the steep slope of the SMF at the massive 
end. 

It is interesting to note that the different cosmologies 
lead also to large systematic errors compared to statistical 
errors. Although many authors in galaxy evolution studies 
claim to account for cosmological parameter uncertainties by 
presenting /i—free results, we recall that, in a flat Universe, 
both flm and Hq enter in the computation of the comoving 
volume and luminosity distances and, in the precision era 
of WMAP and Planck, happen to contribute equally to the 
distance uncertainties. Comparing our results to the recent 
literature is therefore not as simple as scaling the different 
quantities with respect to h, and we must properly account 
for the more complex dependence of distances on flm and 
Ho. 

In comparison, the projected and real-space galaxy clus¬ 
tering as well as galaxy-galaxy leasing are relatively less 
prone to systematic errors. For the effect of cosmology, the 
measurement of projected clustering has no dependence on 
galaxy distances, and the only difference originates from 
the modified galaxy selection caused by the stellar mass 
shift. Interestingly, although the real-space clustering and 
the galaxy-galaxy lensing do depend on galaxy distance mea¬ 
surements, the change in cosmology also has little impact 


at the level of our statistical errors. We can draw similar 
conclusions on the effects of dust extinction modelling and 
photometric calibration. 

Obviously, the stellar mass function is the measured 
quantity suffering from the largest systematic error contri¬ 
bution, compared to the statistical errors. In particular, we 
will see in Section [5] that most of the constraints on the cen¬ 
tral galaxy M* — Mh relationship emanate from the stellar 
mass function and taking into account these systematic un¬ 
certainties when comparing our results with the literature is 
necessary. 

Ideally one would like to estimate a best-fit model 
for each of the re-computed quantities. Unfortunately this 
would be computationally very expensive. Instead, we cre¬ 
ate two sets of measurements: a “statistical error” set based 
on our jackknife error estimate and a “total error” set for 
which we add in quadrature the systematic errors (assuming 
they are Gaussian distributed) and the statistical errors. We 
present in Section [5] separate results for both. 

3.6 Impact of photometric redshift uncertainties 

The dispersion of photometric redshifts may also cause sys¬ 
tematic effects of several kinds, firstly on the stellar mass 
function, as a contribution to the stellar mass scatter, which 
shifts towards higher masses the high-mass end where the 
slope is steep, an effect known as Eddington bias. Secondly, 
the projected clustering amplitude is biased low due to the 
scattering of galaxies falling outside the mass bins. 

We will see in Section 2] that our model properly ac¬ 
counts for these systematic effects caused by photometric 
redshift dispersion, through the parameterisation of the stel¬ 
lar mass scatter. However, catastrophic failures and photo¬ 
metric redshift biases may be more problematic. We have 
demonstrated in Section [2]5] that our catastrophic error rate 
was not h igher than 4% , and based on results from Sec¬ 
tion 3.2 of ICoupon et al.l (l2012| j. such a low contamination 
rate should have no impact on clustering results at our 
statistical error level. To check this statement on the cal¬ 
ibration sample (which means the conclusions are limited 
to the photometric sample with similar properties to the 
spectroscopic sample), we use the VIPERS galaxies with 
spectroscopic redshift and re-compute all stellar masses, as 
well as each observable, using the corresponding photometric 
redshift. We show the measurements in Fig. [S] (solid lines) 
divided by the reference measurement made with spectro¬ 
scopic redshifts and where the error bars are from the sta¬ 
tistical jackknife estimator. From left to right we display the 
results for the stellar mass function, the projected clustering, 
and the galaxy-galaxy lensing signal, all in the mass range 
10^° < M*/M0 < 10*^^ and redshift range 0.5 < z < 1. 

We conclude that for galaxies with similar properties 
to VIPERS galaxies, none of the observables measured with 
photometric redshifts display a large bias with respect to the 
spectroscopic redshift ones. This represents a reassuring san¬ 
ity check for the calibration procedure. Only the projected 
clustering presents a slightly low systematic value, expected 
from the dispersion of redshifts and accounted for in the 
model, through the projection of the modelled 3D cluster¬ 
ing on the redshift distribution constructed from the sum of 
photometric redshift PDFs (assuming that estimated PDFs 
are representative of the true PDFs). 
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Figure 4. Systematic errors affecting the galaxy stellar mass function (top left), the projected correlation function (top right), the 
real-space correlation function (bottom left), and the galaxy-galaxy lensing signal (bottom right). In each panel the gray area symbolises 
the envelope (maximum value) of the re-computed measurement compared to the reference. The error bars are statistical errors from 
the internal jackknife estimator. The “Total” panel represents the symmetric sum in quadrature of all three contributions. Here we only 
show the most massive bins for the clustering and lensing measurements, however we repeated the tests in all mass bins. 
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Figure 5. Measurements made with photometric redshifts divided by those made with spectroscopic redshifts. From left to right: the 
stellar mass function, the projected clustering and the galaxy-galaxy lensing signal, all made with VIPERS galaxies in the mass range 
10^^ < M^/Mq < 10^^ and redshift range 0.5 < z < 1. Error bars represent the statistical error estimates from jackknife resampling. 
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4 MODEL AND FITTING PROCEDURE 

We use the HOD formalism to connect galaxy properties to 
dark matter halo masses. Here we assume that the number 
of galaxies per halo is solely a function of halo mass, split 
into central and satellite contributions. The fitting proce¬ 
dure then consists of finding a set of parameters to describe 
the HOD that best reproduces the observables. 

A key ingredient of the HOD model is the statistical 
description of the spatial distribution of dark matter. We 
assume that the matter power spectrum, the halo mass func¬ 
tion, and the dark matter halo profile are all known quanti¬ 
ties over the scales and redshift range (0.5 < 2 < 1) explored 
in this study. All the technic al details abou t the h alo model 
are given in Appendix A.l of ICoupon et al.l (l2012l i. with the 
exception of the large scale halo bias, f or which we use in th is 
study the fitting formula proposed bv I Tinker et al.l ll2010li . 

The exact way to parameterise the HOD is often at 
the origin of the differences betw een HOD studies in the 
literature. In this paper we follow iLeauthaud et all il201ll i 
who adopted two advanced features: 

• the HOD is a conditional function of the stellar mass 
given the halo mass (this formalism is an extension of 
the conditional luminosity functio n technique developed by 
lYang. Mo fc van den Boschl l2003l l. In this formalism, the 
central galaxy M* — Mh relationship is a parameterised func¬ 
tion representing the mean stellar mass given its host halo 
mass, (M*|Mh); 

• all observables, namely the stellar mass function, the 
projected clustering, the real-space clustering, and the 
galaxy-galaxy lensing signal are fitted jointly. 


4.1 The stellar-to-halo mass relationship 

To describe the central galaxy M* — Mh relationship, we 
adopt the parameterised function /sm-hm proposed by 
iBehroozi, Conroy fc Wechslerl (l2010ll , and defined via its in¬ 
verse: 

logio(/sM-HM) — loglo('^h(M*)) = (7) 



Ml controls the scaling of the relation along the halo mass 
coordinate, whereas M*_o controls the stellar mass scaling. 
/?, S and 7 control the low-mass, high-mass, and curvature 
of the relation, respectively. 


U'logM*(Af*) — UlogM*,0 


f Ml 




(9) 


to account for the change in intrinsic stellar mass dispersion 
as a function of stellar mass. 


4.3 The satellite occupation function 

We describe the satellite HOD for a threshold sample M* 
with a simple power law as a function of halo mass Mh: 

(Wat(Mh|M:)> = 

for which we fix the cut-off mass scale such that 
Mcut = /sm-hm('^*) ■ (11) 


This assumpt i on is based upon the values reported by 
ICouDon et al.l (l2012f) for their equivalent parameter “Mq”. 
We have checked that the exact parameterisation of Mcut 
had very little importance compared to the other parame¬ 
ters and did not change any of our con clusions, in agreemen t 
with the l oose constraints observed bv ICoupon et al.l ll2012ll . 

As in ILeauthaud et al.l (l201l[l . the normalisation Msat 
of the satellite HOD follows the halo mass scaling driven 
by the central M* — Mh relationship, with some degree of 
freedom controlled by a power law: 


Msat _ „ / /sM-Hm(-^*) a 

1012 Mo ^ 1012 Mo J 


( 12 ) 


4.4 Total occupation functions and observables 

Finally, the total HOD is 
(Atot(Mh|M(E)) = 

(Ace„(Mh|M:)) -h (Wat(Mh|M:)), (13) 

and since our measurements are made in bins of stellar mass, 
we transform the threshold HOD functions into binned func¬ 
tions by writing: 

{Ntot{Mi,\Ml\Ml^)) = 

(Atot(Mh|M:i)) - (iVtot(Mh|M:2) . (14) 

Equivalent relations hold for central and satellite binned 
HODs. 

The stellar mass function, the projected two-point cor¬ 
relation function, the real-space correlation function, and 
the galaxy-galaxy lensing signals are computed from the halo 
model and the HOD as detailed in Appendix m 


4.2 The central occupation function 

For central galaxies contained in a threshold sample (M* > 
Ml), the HOD is defined as a monotonic function increasing 
from 0 to 1, with a smooth transition centered on the halo 
mass value corresponding to Ml : 

(Acen(MhlM))) = (8) 

1 [l _ erf ( ~ logio(/sM-HM(Mh)) \ 

2 _ \ %/2(TlogM* (M*) J 

The parameter uiog m* , expresses the scatter in stellar mass 
at fixed halo mass, which we parameterise as: 


4.5 Systematic errors in the model 

As detailed in the previous sections, the HOD formalism re¬ 
lies on an accurate description of the dark matter spatial 
distribution. Here we evaluate the impact of our model un¬ 
certainties and assumptions on the best-fit HOD parameters. 
Ideally, one would like to repeat the fitting procedure to test 
each of the different assumptions of the model, but to avoid 
such a time-consuming exercise, we take the simple approach 
of modifying one feature at a time, and tuning the HOD pa¬ 
rameters by hand to reproduce the modelled quantities de¬ 
rived from the best-fit parameters reported in Section[5] We 
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Table 3. Estimated systematic errors from the model on the central halo mass, log^gMi, and the satellite normalisation, Bsat- The 
total is the sum in quadrature of the errors.. 


Assumption 

Error on logj^Q 

- 12.7) 

loll.5 

Error on Bsa,t (~ 
M*(M0) = lOio 

10) 

IOII.5 

Affected quantities 

0-8 

0.05 

0.05 

1 

0.5 

SMF, clustering (small &: large scales) 

6(Mh) relation 

0.08 

0.1 

— 

— 

clustering (large scale) 

Assembly bias^ 

< 0.04 

< 0.04 

~ 1.5 

~ 1.5 

SMF, clustering (small Sz large scales) 

c(Mh) relation 

0.11 

0.03 

0.1 

0.4 

clustering (small scale), lensing 

Satellite concentration 

— 

— 

1.1 

0.9 

clustering (small scale), lensing 

Stripped sub-halos 

0.09 

0.07 

— 

— 

lensing 

Total 

0.17 

0.14 

2.1 

1.9 

all 


^ Prom^entnCT, Hearin fc van den Bosch (201^. 


explore two stellar mass bins (M* = and 

we focus on the two parameters Mi and Bsat, controlling 
the halo-mass scaling of the M* — Mh relationship, and the 
normalisation of the satellite HOD, respectively. The results 
are shown in Table |3l and we detail below our calculations 
for each assumption listed. 

The power spectrum normalisation parameter, as, is 
currently known to a precision of a few percent. This pa¬ 
rameter has a strong impact on the large scale galaxy clus¬ 
tering, and a larger value would lead to an increased number 
of massive structures, hence shifting the massive end of the 
halo mass function towards more massive halos. Choosing 
Planck over WMAP7 cosmology (as for the tests in Sec¬ 
tion I3.5II . would result in a 5% increase in ag, leading to 
relatively small changes in best-fit HOD parameters, of or¬ 
der of a few percent. 

Halo bias uncertainties originate from the measurement 
of the bias-to-halo mass relation b{Mh) using n-body sim¬ 
ulations, affected by low-mass resolution, small volume, or 
the limitations of halo identihcation techniques. In the low- 
clustering regime, the typical errors on the bias are as small 
as a few percent (iTinker et al.l 1201(11 , however the rather 
shallow slope of bia s versus halo mass (see e.g. Fig. 18 of 
ICouDon et al.ll2012ll translates into a larger uncertainty in 
the deduced halo mass. In the high-mass regime, errors are 
mainly dominated by the sample variance of simulations, up 
to ~ 10%, but have fewer impact on the deduced halo mass 

owing to the steeper slop e in this regime. _ 

_ The assembly bias ijZentner, Hearin fc van den Boschl 

l20l4 and references therein) refers to the correlation be¬ 
tween clustering amplitude and time of halo formation, 
whereas in our model the bias is assumed to vary only with 
halo mass. The effect is stronger when selecting a popula¬ 
tion of galaxies based on a parameter correlated with halo 
formation history, such as the star formation rate, but mod¬ 
erate when considering the full galaxy population selected 
by stellar mass only. In this case, and in the mass regime 
explor ed in this study, IZentner. Hearin fc van den Boschl 
(l2014ll found that the systematics caused by assembly bias 
on HOD parameters do not exceed 10% to 15%. 

In our model, the dark matter halo pr o hle is 
assumed to follow a iNavarro. Frenk fc White! (|l997l ~l 
(NFW) profile. While leasing observati o ns tend to 
favour NFW prohles llUnigtsu_et-aL l201l|: lOkabe et al.l 
120131 : ICouDon. Broadhurst fc Umetsu l2013ll . the mass- 
concentration relation - driving the slope of the pro¬ 
file - remains uncertain. We have used a simple mass- 


concentrati on relation based on theoretical predictions (up¬ 
dated f romiTakada fc Jain l2003h and empirical redshift evo¬ 
lution dBullock et al.l 2001 ), but more re c ent re lations such 
as the work from Munoz-Cuartas et al.l ll2010h have been 
measured. C ompared to our conce ntration values, the dif¬ 
ference with iMunoz-Cuartas et ^ rises from 11% at Mh ~ 
10^^Mq to 30% at ~ 1O'^®M0 (with a minimum of 2% at 
~ 10^^ Mq). These systematics affect the slope of the small 
scale clustering and galaxy-galaxy leasing. We estimate that 
if all of our constraints came from lensing, this may result 
in a 28% systematic error in Mi. 

We assume that the satellite distribution in the halo fol¬ 
lows the dark matter density profil e. However, this a s sump - 
tion may not be always true and iBudzvnski et al.l (l2012fl 
tested this hypothesis from a stacked analysis of massive 
clusters from the SDSS. They found a typical factor of 2 
(with ~50% scatter) lower concentration of th e satellite dis¬ 
tribut ion compared to dark matter, whereas iMuzzin et al.l 
(l2nnitl measured a value closer to dark matter around 2 ~ 
0.3, and Ivan der Burg et al. I ||2014) a relatively high concen¬ 
tration at 2 = 1. These trends may show a redshift evolution 
of the concentration or can simply be inherent to the diffi¬ 
culty of observationally measuring the satellite distribution. 
In Table [3] we report the impact on Hsat after setting the 
satellite concentration a factor two higher than that of dark 
matter. The effect on Bsat does not exceed 11%. 

Finally, in our model we neglect the lensing contribu¬ 
tion of the sub-halos h osting the satellite galaxie s. This ef¬ 
fect, first introduced bv iMandelbaum et al. (|2005h under the 
term “stripped satellite central profile”, assumes that a frac¬ 
tion of the satellite halos survive inside the host halo and 
further contribute to the lensing signal at small scales. As a 
result, the lensing contribution of those sub-halos adds up to 
the central-galaxy halo term in such a way that the best-ht 
host halo mass gets reduced compared to a model in which 
the c ontribution of sub-halos is neglected. iHudson et al.l 
1I2OI5II quantify the systematic change in best-fit halo mass 
as a systematic decrease by a factor of ~ (1 + /sat), where 
/sat is the fraction of satellites in the sample. Assuming a 
satellite fraction between 20% and 30%, this leads to a sys¬ 
tematic error of up to 0.09 in log^Q Mi. This number must be 
read as if all the constraints would come from lensing only. 
In our study where the stellar mass function and the clus¬ 
tering signal-to-noise ratio is higher than that of the lensing, 
this effect plays relatively little role, and our results would 
not significantly change if we accounted for it. 

The sum in quadrature of these model systematics is 
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shown as “Total” in Table [3] Intermediate-stellar mass bins 
(~ 10 ^°Mq) seem to be most affected, with an error of 0.17 
for logj^g Ml (~ 50% in Mi) and 2.1 (~ 20%) for Bsat. We 
will see below that these values dominate over the typical 
statistical and systematic errors from the measurements in 
this mass regime. However, as each of these systematic er¬ 
rors affects the observables in a different way and we fit all of 
them jointly, one must see these numbers as pessimistic es¬ 
timates. The high-stellar mass bin (~ Mq) is equally 

affected but in a regime where statistical errors are large, 
hence leading to a smaller impact. 


4.6 MCMC sampling 


We write the combined log-likelihood as the sum of each 
observable 

— 2 InZl = xip(rp) + Xae i (15) 

spl spl spl 


where individual are computed as: 


2 

X = 


Et 


Xi-Xi 






(16) 


using the covariance matrices evaluated for each measure¬ 
ment as described in Section[3](X and X represent the mea¬ 
sured and modelled observables, respectively). Each observ¬ 
able is summed over the samples (“spl”) as described in 
Table [2] The “i” and “j” subscripts refer to the stellar mass 
(stellar mass function) or transverse separation (clustering 
and lensing) binning of each measurement. 

We find the best-fit parameters and posterior distri¬ 
bution (assuming flat priors for all parameters) employing 
the Markov Chain Monte Carlo (MCMC) sampling tech¬ 
nique with the Metro polis-Hastings sam pler from the soft¬ 
ware suite Cosmopmc (I Wraith et al.ll20C)^ . We check for in¬ 
dividual chmr]_convergmceand chain-to-chain mixing us¬ 
ing the ICelman fc RubinI ill992l l rule from the R-language 
CODA packagj^°l. We find a typical chain-to-chain mixing 
coefficient (potential scale reduction factor) to be equal to 
1.01, and the acceptance rate around 30%. 

In practice, we first evaluate a diagonal Fisher matrix 
at the ma ximum likelihood point found using the Amoeba 
algorithm dPress et al. IH) and run 10 chains in parallel 
with the inverse Fisher matrix as the MCMC sampler covari¬ 
ance matrix. The acceptance rate is usually very low due to 
the noisy diagonal Fisher matrix affected by some strong 
correlations between parameters. Once the chains have con¬ 
verged (after typically 5 000 — 10 000 steps) we compute the 
final likelihood covariance matrix after rejecting the burn-in 
phase of the chains (a few thousand steps). This covariance 
matrix is used as the input sampler covariance matrix of a 
second and final MCMC run, in which 10 chains of 30 000 
steps each are computed in parallel and combined together 
assuming a burn-in phase of 2 000 steps and checking for 
proper mixing. 

We run the full MCMC procedure twice. The first run 
is performed using the statistical covariance matrices from 
the jackknife estimator and the second MCMC run uses the 


http://cran.r-project.org/web/packages/coda/citation.html 


Table 4. HOD best-fit parameters and 68% confidence limits 
(CL) for the statistical errors (top) and total errors (bottom). 


Jackknife resampling errors 


Parameter 

mean 

upper CL 

lower CL 

logic ^1 

12.84 

0.020 

-0.026 

logic ^*,0 

10.98 

0.015 

-0.019 

p 

0.48 

0.017 

-0.021 

5 

0.63 

0.094 

-0.073 

7 

1.60 

0.166 

-0.202 

'riogM*,0 

0.337 

0.045 

-0.035 

A 

0.21 

0.047 

-0.044 

-^sat 

10.87 

0.443 

-0.416 

/^sat 

0.83 

0.038 

-0.035 

a. 

1.17 

0.020 

-0.021 



Total 

errors 


Parameter 

mean 

upper CL 

lower CL 

logic 

12.67 

0.124 

-0.083 

logic ^*,0 

10.90 

0.082 

-0.067 

p 

0.36 

0.077 

-0.051 

& 

0.75 

0.193 

-0.151 

7 

0.81 

0.477 

-0.386 


0.394 

0.100 

-0.074 

A 

0.25 

0.082 

-0.083 

-^sat 

9.96 

0.938 

-0.845 

/3sat 

0.87 

0.078 

-0.065 

a. 

1.14 

0.040 

-0.038 


total error covariance matrices, which are constructed from 
the statistical covariance matrices plus the systematic error 
estimates added in quadrature to the diagonal, as described 
in Section 13.51 


5 RESULTS 

Best-fit parameters with 68% confidence intervals are given 
in the top panel of Table|3]for the statistical- and total-error 
MCMC runs. The ID and 2D likelihood distributions are 
shown in Fig. EH The reduced xS for the statistical-error 
fit is X^/(-^points - iVparameters) = 260/(160 - 10) = 1.7, 
which is an overestimate given the correlations neglected 
in the computation of the log-likelihood. Firstly, we recall 
that the lensing and clustering measurements are affected 
by a sample-to-sample correlation due to the scatter in stel¬ 
lar mass. The re-use of background galaxies in the lensing 
measurements causes an additional sample-to-sample cor¬ 
relation. Secondly, the projected and real-space clustering 
are correlated, as both observables bring similar informa¬ 
tion. This mostly affects the satellite distribution parameter 
errors, which could be slightly underestimated. Finally, the 
few number of sub-samples (64) used in the computation 
of a noisy covariance matrix may have biased the inverse 
estimate and contributed to an increase in xS- 


5.1 Measurements and best-fit models 

The measured stellar mass function and best-fit model are 
displayed in Fig. (6] Statistical error bars and correspond- 
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Figure 6. Measured stellar mass function and best-fit model in 
the range 0.5 < 2 : < 1. The statistical errors from the jack¬ 
knife estimate are shown as black thick lines, whereas the to¬ 
tal (statist ical plus systema tic) error bars as dotted lines. The 
COSMOS Gilbert et al.l20l'3) and VIPERS llDavidzon et al.l2Q13tl 
mass functions are displayed with their respective statistical er¬ 
rors as shaded areas. 

ing best-fit model are shown as thick black lines, whereas 
total (statistical plus systematic) errors and corresponding 
best-fit model are represented in dotted lines. We compare 
our measurements with the COSMOS mass function eval- 
uated in th e rang es 0.5 < 2: < 0.8 and 0.8 < 2; < 1.1 by 
lllbert et al.l (|201^. an d the VIPERS stellar mass function 
(iDavidzon et al.luOlSt) . measured in the range 0.5 < z < 1 
(1. Davidzon, private communication). 

The clustering measurements and best-fit models are 
shown in Fig. [T] The projected two-point correlation func¬ 
tions w{6) are displayed in the top panels. The mass 
ranges are given in each top-right panel corner in units of 
log(M*/M0). Similarly, the real-space two-point correlation 
functions W(J'p) are displayed in the bottom panels. 

The galaxy-galaxy lensing measurements and best-fit 
models are shown in Fig. [S] The most massive lensing bin 
features a few data points lower than the model around the 
transition between the central and the satellites term. 

For all observables, we report good agreement between 
the data and the model. The constraints on the shape of the 
central M* — Mh relationship (parameterised by logj^o-^i) 
logj^Q Mstaro, /?, 7 and 5), are mostly driven by the high 
signal-to-noise stellar mass function measurements. Satellite 
HOD parameters (Bsat, /3sat and q) are mainly constrained 
by the clustering and lensing measurements. The amplitude 
of clustering at small scale is directly proportional to the rel¬ 
ative number of satellites, hence giving strong leverage on 
the satellite galaxy HOD. Additional information is given 
on scales r ~ 0.1 Mpc from lensing, through the satellite 
lensing signal. The dispersion in M^, at fixed M^, parame¬ 
terised in amplitude by criogM. ,o and in power-law slope by 


A, is mainly constrained by the high-mass end of the stellar 
mass function and the amplitude of the galaxy-galaxy lens¬ 
ing signal in the most massive bins, resulting in a high-mass 
(M* ~ 10 ^^Mq) scatter of approximately niogM* — 0.2 in 
both the jackknife and total error cases, and a medium mass 
(M* ~ 1O^°M0) scatter of criogM* — 0.35. 

Because the stellar mass function is most affected by 
the inclusion of systematics in the error budget, we note a 
significant increase in uncertainties associated with the pa¬ 
rameters driving the central M* — relationship. From 
Table |4l we report an increase from a factor ~ 3 in the er¬ 
ror in 7, up to a factor ~ 6 in the error in logjQ Mi. HOD 
parameters describing the satellite occupation function such 
as Bsat , P OT a show substantially less sensitivity to the ad¬ 
dition of systematic errors in the error budget (a maximum 
of factor ~ 2 increase is found). This is explained by the 
relatively smaller contribution of systematic versus statisti¬ 
cal errors affecting the clustering and lensing measurements, 
compared to the stellar mass function. 

The occasional large differences between best-fit param¬ 
eters from statistical alone and total errors, seen in Table IH 
do not lead to significantly different derived quantities, ow¬ 
ing to the strong correlations between parameters. This is 
confirmed by the almost indistinguishable dotted lines and 
thick lines in Figs[EHll and is most probably a consequence 
of having symmetrically added the systematic errors to the 
statistical errors. 


5.2 Central M* — relationship and the SHMR 

In Fig. [9] we show the best-fit central galaxy M* — Mh re¬ 
lationship (left panel) as parameterised by Eq. ([§1), and the 
stellar-to-halo mass ratio (SHMR, right panel). The SHMR 
is shown as function of host halo mass and is derived for the 
central galaxy in dark grey (from the M* — Mh relationship), 
the satellites in light grey (integrated over the galaxies above 
a mass threshold of M* > 10^'^ M0), and the total in black. 

The shaded areas represent the 68% confidence limits, 
and in the bottom left panel we have shown the results ob¬ 
tained with statistical errors in light blue and with total 
errors in black. As for the stellar mass function, the statisti¬ 
cal uncertainties grow by a factor ~ 2 — 4 in the lower mass 
regime, when incorporating systematics. 

The central SHMR peak position is indicated by a black 
arrow located at Mh,peak = 1.92jlQ ^4 x 1O^^M0. The SHMR 
peak value is SHMRpeak = x 10~^. When account¬ 

ing for satellites, the peak position and value do not sig¬ 
nificantly differ from the estimates for centrals only. How¬ 
ever, a remarkable result highlighted in this figure is the 
increasing contribution of stellar mass enclosed in satellites 
as function of halo mass. When reaching cluster-size halos, 
this contribution reaches over 90% (and presumably higher 
when accounting for satellite galaxies with masses lower 
than 1O^°M0). However we stress that we do not take into 
account the intra-cluster light (ICL), which is challenging to 
quantify using ground-based photometric data. 

5.3 Comparison with the literature 

In Figs [To] and [ 11 ] we compare our best-fit M* — Mh relation¬ 
ship for central galaxies with a number of results from the 
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Figure 7. Galaxy clustering measurements (data points with error bars) and best-fit models (thick lines). The top panels show the 
projected w{9) from the photometric sample (the measurements are corrected for the integral constraint), and the bottom panels show 
the spectroscopic real-space Wp(rp). The thick error bars associated with thick lines represent the statistical errors and subsequent best-fit 
models, whereas dotted lines are for total errors. The mass ranges in the top right corner of each panel are given in /Mq). 


literature. As described in Section lU our relation describes 
the mean stellar mass at fixed halo mass which is, due to the 
scatter in stellar mass, not equivalent to the mean halo mass 
at fixed stellar mass. This issue becomes particularly impor¬ 
tant when the slope of the stellar or halo mass distribution 
is steep (i.e. at high mass). Therefore, we have re-computed 
our results using the latter definition and we consistently 
compare our results with the literature in each case. 


When required, we convert ha lo masses to ou r viria l 
definition using the recipe given bv iHu fc Kravts^ (l2003l ) 
in their Appendix C and, following lllbert et al. 1 20ld) . we 
divi de stellar masse s by a factor 1.74 and 1.23 to convert 
from ISal peteil ll 1 9551) and “Diet” Salpeter IMFs, respectively, 


to our 


Chabrierl I MF stellar masses. We apply no correction 


to iKroupX i 200 ll ) IMF stellar masses. 


The mean redshift, measured from the sum of the pho¬ 
tometric redshift PDFs, is found to be {z) = 0.82 for our 


measurements in the range 0.5 < z < 1.0 (M* > Mq) 

and (z) = 0.65 in the range 0.5 < z < 0.7 (10^° < 
M*/M0 < 10^°'"^°). We point out that the lensing signal is 
more sensitive to lower-redshift lens galaxies characterised 
by a higher signal-to-noise (due to the more numerous back¬ 
ground sources), and is likely to be more representative of 
a lower redshift population, but this effect is assumed to be 
small compared to the lensing statistical errors. 


5.3.1 (M*|Mh) results 

We first compare the results for (M*|Mh) in Fig. 1101 The 
black shaded area shows our results for the central galaxy 
relationship with 68% confidence limits from the total errors. 
The total errors consist of the statistical uncertainties from 
jackknife resampling in addition to three sources of system¬ 
atic effects from the measurements: the cosmology chosen 
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Figure 8 . Galaxy-galaxy lensing signal measurements (data points with error bars) and best-fit models (thick lines). As in Fig. 0 thick 
and dotted lines are for statistical and total error results, respectively. The mass ranges in the top right corner of each panel are given 
in log(M*/M 0 ). 



10l2 iq13 ^qU ^q15 

M. [Mq] 



Mh [Mo] 


Figure 9. Best-fit M* — Mh relationship (left) and stellar-to-halo-mass ratio (SHMR, right). The black shaded areas represent the 
confidence limits from the total errors. The bottom left panel shows the confidence limits interval as a function of halo mass in the case 
of statistical errors (from jackknife resampling in light blue) and total errors (in black). The SHMR is derived as function of host halo 
mass for the central galaxy (dark grey), the satellites (light grey) and the sum of both (black). The peak value of the central SHMR is 
indicated by the black arrow. 
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Figure 10. The best-fit M* — Mh relationship for central galaxies, shown in the black shaded area (total-error based 68% confidence 
limits), compared with a number of results from the literature at similar redshifts. The results shown here represent the mean stellar 
mass at fixed halo mass or halo-mass proxy (X-ray temperature or satellite kinematics), (M*|Mh), but plotted Mh as function of M* to 
ease the comparison with the literature. We perform appropriate halo mass conversions and IMF stellar mass corrections when required. 
The length of t he grey arrow in t h e bot tom r ight corner shows the shift ('^ 0.2 dex) meas ured from the direc t comparison between stellar 
masses used in lLeauthaud 3 1 I 2 OI 2 II and iGeorge et al.l ll201lh . compared to those in lllbert et all ll2010h which were estimated in a 
similar way to this study. The error bar on the bottom-right corner indicates the typical systematic uncertainty arising from the model. 


among widely-used ACDM models, the fine-tuning of our 
dust extinction law modelling, and potential biases in the 
photometry/calibration. We recall that this list of system¬ 
atic uncertainties is not exhaustive and, for example, ignores 
the choice of SPS models, which may be responsible for even 
larger systematic effects. An estimate of the systematic er¬ 
rors from the model, as detailed in Section|431 is also shown 
in the bottom-right corner. 

iBehroozi et al.l ll2013l '). shown as the light-blue shaded 
area, put constraints on the M^, — relationship by popu¬ 
lating dark matter halos in simulations and comparing abun¬ 
dances using observed stellar mass functions from a number 
of surveys. They characterised the uncertainties affecting 
stellar mass estimates by accounting for a number of system¬ 
atic errors. In particular, unlike in our systematic errors, the 
authors had to include uncertainties arising from the choice 
of the IMF and the SPS galaxy templates, necessary when 
combining the stellar mass functions from several works us¬ 
ing different stellar mass measurement methods. Here we 
consider their results at z ~ 1. A significant difference with 
our model resides in the assumption that satellite galaxies 
in larger halos are seen as central galaxies in sub-halos. To 
circumvent the difficulty of accurately predicting a sub-halo 


mass function (e.g. complications from tidal stripping), the 
galaxies in sub-halos at the time of interest are matched to 
their progenitors at the time of merging onto the central 
galaxy halo, under the assumption that the M* — Mh evolu¬ 
tion at a given stellar mass is identical whether the host halo 
is isolated or inside a larger halo. In comparison, our model 
is a “snapshot” of the galaxy halo occupation at a given 
time, where the satellite distribution is mainly constrained 
by galaxy clustering. 


The results from iLeauthaud et al.l (l2012l i in COSMOS 
are shown in brown and green at redshifts 2 ~ 0.6 and 
2 ~ 0.9, respectively. We observe a small discrepancy which, 
compared to our results, is unlikely to be explained by dif¬ 
ferences in the modelling of the HOD (since the model is 
essentially identical), nor the sample variance as confidence 
limits do not overlap. A difference in stellar mass estimates 
on the other hand is more likely to be at the origin of the 
discrepancy. To check this hyp othesis, we have co mpared 
the stellar mass estimates from lllbert et all ll2010h . which 
were m easured in a sim il ar wa y to this study, with those 
us ed in Leauthand^tah liioli) with the method described 
in iBundv et al.l 1 2006t l. We measured an offset of ~ 0.2 dex, 
illustrated in Fig. [T0]as the grey arrow. Part of the difference 
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seems to be explained by the separate choice for the dust ex¬ 
tinction law made in each study (which may typically cause 
a ~ 0.14 dex offset, see Section f2.611 . However we note that 
in both cases the same IMF and set of SPS models were 
used, which leaves us without a complete understanding of 
the difference. _ _ 

The results by IWang k, Jind ll2010lf are shown as the 
blue short-dashed line. Their model is based on a HOD mod¬ 
elling of the stellar mass function and rea l-space galaxy clus¬ 
tering where, as in iBehroozi et af] (l2013l f . the treatment for 
satellites is not based on the distribution of sub-halos in the 
host halo but on the M* — Mh relationship at the time of 
infall. 

iMoster. Naab fc Whitel (l2013f l. shown as the red dot- 
dashed line, also used abundance matching and provided a 
redshift-dependent parameterisation of the central M* — Mh 
relationship that we have calculated at z = 0.8. As above, 
the satellites are matched to their halos at the epoch of 
merging. Their relation is in good agreement with ours at 
intermediate mass, however, it shows a steeper dependence 
on stellar masses at higher mass. 

The green dots with error bars ar e from the HOD mod¬ 
elling results of IZheng. Coil fc Zehavil (l2007l l. based on real- 
space clustering and number density measurements. Here 
we show their results for DEEP2, a deep spectroscopic sur¬ 
vey with high density z = 1 galaxies. Without deep NIR 
data, the authors have computed mean approximate stel¬ 
lar masses for galaxy samples selected in bins of luminosity. 
This source of uncertainty is not shown on the plot, however 
one may expect a large scatter and potential biases due to 
this conversion. 

The o range bow-t i es wit h error bars represent the re- 
sultfEB by IWake et all (1201 ill in the NEWFIRM Medium 
Band Survey at redshift 2 ~ 1.1, from the combination of 
NIR-selected galaxy clustering and number density measure¬ 
ments. Their results are in good agreement with ours. 

The five next results were produced using galaxy clus¬ 
ter sam ples associated with their brightest cluster galaxies 
fBCGl. ICeorge et al.l (l201l[l built up a catalogue of central 
versus satellite galaxies in COSMOS, matched to an X-ray 
detected grou p/cluster sample with r obust halo masses from 
weak lensing (iLeauthaud et aDl2009o . From their catalogue 
we have computed the mean of stellar mass and halo mass 
values for clusters in the range 0.5 < 2 < 1, shown as the 
single red triangle (the error bars show the standard de¬ 
viation in halo a nd stellar masses ) . As t hey used identical 
stellar masses to ILeauthaud et all (I 2 OI 2 I I , we also expect a 
systematic difference in stellar masses compared with our 
estimates. 

From Sunyaev-Zel’dovich (SZ) detected clusters using 
the Atacama Cosmology Telescope, iHilton et al.l 1 I 2 OI 3 II pre¬ 
sented the measurements of the galaxy properties between 
0.27 < 2 < 1.07. Member galaxies were identified from high 
density spectroscopic observations, and stellar masses were 
measured from Spitzer IRACl-2 mid-infrared (MIR) fluxes. 
Halo masses were estimated from satellite kinematics. Here 
we show the mean halo mass versus mean BCC stellar mass, 


Here we use updat ed results compared to the original publi¬ 
cation, estimated with iBruzual fc Charlotl (l20n3l i templates and 
with rectified /i-scaling (D. Wake, private communication) 


represented by the single blue dot with errors bars (stan¬ 
dard deviations of both masses). Their results appear to be 
in good agreement with our M* — Mh relationship, although 
our constraints on such high-mass clusters are extrapolated 
from the few clusters more massive than 4 — 5 x IO^^’Mq 
expected in our sample. 

We show as a single light blue diamond the 
mean halo mass v ersus mean BCC stellar mass from 
Ivan der Burg et al.l (l2014l ) in the CCLASS/SpARCS clus¬ 
ter sample at 2 ~ 1. Calaxy cluster members were identified 
from intensive spectroscopic observations, and halo masses 
were estimated from satellite kinematics. We note that stel¬ 
lar masses were measured from a similar combination of 
data, redshift range and volume size as ours, however the 
methodology used to link halo mass to galaxy stellar masses 
was rather different. Thus, the agreement with our high- 
mass M* — Mh relationship within the sample variance is 
quite remarkable. 

Results from iBalogh et al.l ll2014h are shown as the 
downward purple triangles. Halo mass measurements were 
made using satellite kinematics for a sample of 11 
groups/clusters in the COSMOS field. We show the mean 
and standard deviation of their measurements split into two 
halo mass bins (the 11 groups are split into 5 and 6 groups 
below and above Mh = 9x 10^^Mq, respectively). Although 
their results suffer from large sample variance, they are in 
broad agreement with our results and with the rest of the 
literature. 

Finally, the single red square with error bars shows 
the mean of halo mass measurements from a weak lens¬ 
ing analjrsis_ofX-ray_selected clusters in the CFHTLenS 
by iKettula et al.l |20M^ versus the mean st ellar mass of 
associated BCGs (iMirkazemi M. et al.llsubm.h . We have re¬ 
measured stellar masses of those BCGs in a consistent way 
to this study (with the exception of missing NIR data for 
most of the BCGs, which may increase the scatter in stellar 
mass). Despite the lower redshift range, the identical pho¬ 
tometry and lensing catalogue makes the comparison rel¬ 
evant to our results, where the expected difference should 
arise solely from redshift evolution, although the large sta¬ 
tistical uncertainties prevent us from drawing strong conclu- 


5.3.2 (Mh|M*) results 

We compare the results for (Mh|M*) in Fig. 1111 To express 
the mean halo mass at fixed stellar mass (Mh|M*) from our 
results, we derive it from the mean stellar mass at fixed halo 
mass (M*|Mh) using the Bayes theorem relating conditional 
probability distributions: 


P{Mh\M0 oc P(M*|Mh) X P{Mh). 


(17) 


We can then compute (Mh|M*) as the expectation value of 

P(Mh|M*): 


(Mh|M*) 


/ P{M.\Mh)P{Mh)MhdMh 

f P(M*|Mh)P(Mh)dMh 


(18) 


with 


P(M*|Mh) 


d(Ncen(Mh|M„)) 

dM* 


(19) 


the distribution of central galaxies given a halo mass, and 
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Figure 11. The best-fit Mir — Mh relationship for central galaxies, shown in the black shaded area (total-error based 68% confidence 
limits), compared with a number of results from the literature at similar redshifts. Unlike in Fig. 1101 the results shown here represent 
the mean halo mass at fixed stellar mass (Mh|M*). W e perform appropriat e halo mass conversions and IMF stellar mass corrections 
when required. The relatively low halo masses found by [Hudson et all ll2015l) is linked to a different treatment of the satellite sub-halo 
contribution to the lensing signal at small scale (see text for details). 

= ( 20 ) 

the halo mass function. _ _ 

We show the results of iFoucaud et all (|201fll l at 2 ~ 1 
from clustering measurements in the UKIDSS-UDS field as 
the blue squares with error bars. The UKIDSS-UDS held 
is a small patch of ~ 1 deg^ with deep NIR and optical 
data. They have converted their clustering amplitude mea¬ 
sured in bins of stellar mass into halo masses, using the ana¬ 
lytica l galaxy-bias halo-mass relationship from I Mo fc Whita 
( 19961) . As they do not use any constraints from galaxy num¬ 
ber density, their error bars are dominated by sample vari¬ 
ance and uncertainties on the projected galaxy clustering. 

Green upward triangles represent the results by 
IConrov et al.l (120071 '). Halo masses were derived from satel¬ 
lite kinematics using spectroscopic measurements from the 
DEEP2 survey. Since the authors have selected their samples 
based on bins of stellar masses, we can compare their results 
with our M* — Mh relationship. The agreement is 

found to be good. 

Results from clusteri ng measurem e nts in the CEHTLS- 
DEEP/WIRDS helds by iBielbv et al.l ll20l4 ) are displayed 
by the brown bow-ties with error bars. We select all mass bin 
results in the range 0.5 < 2 < 1. Although the total field-of- 
view is small (~ 2.4 deg^), the combination of four indepen- 


dent fields allow ed them to reduce the cosmic variance. As in 
iFoucaud et al.l . they used an analytical prescription based on 
the large-scale clustering amplitude to estimate halo masses 
per bin of stellar mass, so that their results should be com¬ 
pared to our (Mh|M*) results. The two points well above 
the other results correspond to the measurements at 2 ~ 0.7 
and seem to disagree with our constraints and the rest of 
the literature. The authors claim to have observed an un¬ 
usually high clustering signal at those redshifts, potentially 
explained by cosmic variance effects. 

Results by iHevmans et al.l (l2006l l in the COMBO- 
17/GEMS field are shown as the downward light-blue tri¬ 
angle with error bars. Here we have picked their unique 
measurement at 2 > 0.5. Halo masses were measured us¬ 
ing weak lensing with galaxy shapes from the Hubble Space 
Telescope observations. 

We sho w as red diamon d s the results for 2 ~ 0.5 red 
galaxies by Ivan Uitert et al.l (l201lh in the Red Sequence 
Gluster Survey 2, a medium-deep GFHT-MegaGam survey 
in three bands (gri) which overlaps 300 deg^ of the SDSS. 
The authors have measured the galaxy-galaxy lensing signal 
for SDSS lens galaxies with a spectroscopic redshift using 
background source galaxies from the RGS2 survey. Here the 
large area permits a high signal-to-noise measurement for 
very massive galaxies from lensing only. Their results are 






▲ 

M 

T 

♦ 

• 

A 

O 

A 


CFHTLS-NIR/VIPERS, z~0.8 (this study) 

Foucoud et ol. (2010), z~1.0 
Conroy et ol. (2007), z~0.8 
Bielby et ol. (2013), z~0.8 
Heymons et ol. (2006), z~0.7 
von Uitert et ol. (2011), red, z~0.5 
Velonder et ol. (2014), red, z~0.3 
Velonder et ol. (2014), blue, z~0.3 
Hudson et ol. (2015), red, z~0.7 
Hudson et ol. (2015), blue, z~0.7 


i 







Typical model 
systemotics 


© 2014 RAS, MNRAS OOO.ITltgTI 









































HOD in CFHTLenS/VIPERS 21 


consistent with ours as this mass bin (> 3 x lO^^M©) is 
dominated by red galaxies. _ 

W e compare our results with those from IVelander et al.l 
(l2014ll at z ~ 0.3, shown as filled symbols (red dots and 
blue trian gles for red and blue galaxies, respectively), and 
those from [Hudson et ^ ll2015l i at 2 ~ 0.7, shown as empty 
symbols (red dots and blue triangles for red and blue galax¬ 
ies, respectively). In both studies, halo mass measurements 
were obtained from galaxy-galaxy leasing measured using 
the CFHTLenS leasing catalogue and stellar masses com¬ 
puted in a similar way to this study, with the exception that, 
in both cases, no NIR data were available at the time. Th is 
mostly affects the ste llar mass estima tes of iHudson et al.l at 
2 ~ 0.7 which, unlike IVelander et al.l at 2 ~ 0.3, do not ben¬ 
efit from the leverage of the CFHTLS 2 -band. We expect the 
M* — Mh relationship of the full galaxy population to lie be¬ 
tween tho se of the red an d blue populations, however the re¬ 
sults from iHudson et ^ lie below our results for both galaxy 
populations. The bias caused by the scatter in stellar mass 
partially explains this difference (by shifting t heir mean stel- 
lar mass to higher values), but not entirely: [Hudson et al.l 
account for the contribution of sub-halos around satellites 
occurring at small scale in the leasing sign al, whereas we do 
not (see Section ri.dPh . As IVelander et all also accounted for 
sub-halos in their leasing model, we cannot exclude that the 
apparent good agreement may result from a redshift evo¬ 
lution going in the opposite direction, and requires further 
investigation. 

5.3.3 The total SHMR 

In Fig. [ 12 ] we show the SHMR as function of halo mass 
compared with observations from the literature. The black 
shaded area represents the total SHMR as the sum of the 
central and satellite contributions. The central SHMR (in 
dashed line on the figure) is simply derived from the central 
M* — Mh relationship. The satellite SHMR (in dot-dashed 
line on the figure) is computed from the sum of satellite stel¬ 
lar masses over the halo occupation function at each halo 
mass, with a lower integration limit of M* = 10^°M©. The 
total baryon fraction compared to dark matter in the Uni¬ 
verse is assumed to be 0.171 and r epresented on the fi gure 
by the grey shaded area on the top dPunklev et al.l[2009l . the 
width of the line represents the uncertainty). 

In green we display the total SHMR from 

iLeauthaud et all (I 2 OI 2 I I measured at 2 ~ 0.9. The 

procedure to compute the total SHMR is identical to ours, 
i.e. the integrated stellar masses from the satellite HOD 
were added to the central stellar mass at each halo mass. 
The authors adopted a mass threshold of 10®'®M©, which 
does not change the integrated stellar mass from satellites 
by a large amount compared with a cut of > 10^® M©. As 
shown in Fig. 1101 part of the vertical shift is explained by 
the systematic difference in stellar mass estimates. 

We show in light blue the central SHMR from 
Behroozi, Wechsler fc ConrovI (l2013l 'l. As seen in Fig. llOl the 
agreement with our central SHMR is good, although their 
peak is located at a slightly lower halo mass value. 

This point is al so investigated in detail in Appendix D of 
[Hudson et al.l l[2015l l. 


_ The red triangle shows the results by [George et al.l 

|[201llf in COSMOS in the redshift range 0.5 < 2 < 1. 
The point represents the mean total stellar mass divided 
by the halo mass versus the halo mass, and the error bars 
the standard deviation in each direction. Here we computed 
the total stellar mass as the sum of the central galaxy stel¬ 
lar mass plus the stellar masses of associated group mem- 
bers with Mi, > 10^®. As they used the stellar masses of 
[Leauthaud et all , the agreement is consistently good with 
their results, however shifted compared to ours. 

The single blue dot with error b ars marks th e mea n 
and standard deviation of estimates bv iHilton et al.l l[2013ll . 
Here the total cluster stellar mass is measured from the 
background-subtracted sum of galaxy IRAC fluxes within 
f ?500 from the BCG. Based on the stellar mass completeness 
computed by [llbert et al.l l[2010[ ). a IRAC AB magnitude 
cut of 24 gives a complete passive galaxy sample down to 
M* = 10®M© at 2 ~ 0.5. With an IRAC completeness AB 
magnitude li mit of 22.6, it is therefore safe to assume that 
[Hilton et al.l are complete above 10^®M© at 2 ~ 0.5, which 
matches our sample. We then conclude that their measure¬ 
ments are in good agreement with our result s. 

Results from [van der Burg et al.l ([2014h are shown as 
the single light-blue diamond, representing the mean SHMR 
versus halo mass with its standard deviation. Total stellar 
masses are computed as the sum of the BCG stellar mass 
and the stellar mass from galaxy members spectroscopically 
identified and corrected for target sampling rate. The au¬ 
thors have checked that for > 10^® M© galaxies, which con¬ 
tribute the most to the total SHMR (see their Fig. 2), the 
spectroscopic success rate reaches 90%. We note that the 
median stellar mass completeness ~ 10^®’^®M© is slightly 
higher than ours (limited by their Jfs-band data), however 
the contribution of satellites compared to a mass limit of 
10^® M© will not significantly change the total SHMR and 
their measurements can be fairly compared to our results, 
and we observe an excellent agreement. Interestingly, the au¬ 
thors conclude that when comparing with the literature, no 
redshift evolution in the total SHMR at high mass is found 
below 2 ~ 1 a nd the compa rison with our results (2 ~ 0.8) 
and those from iHilton et aH (2 ~ 0.5) confirm their findings. 

The two purple downw ard triangles represent the results 
from iBalogh et al.l (120141 ') in the GEEC2 survey in COS¬ 
MOS. Here we show the mean and standard deviation of 
the SHMR versus halo mass in two halo mass bins. Galaxy 
members are identified from the spectroscopic redshift when 
available or using the PDF-weighted photometric redshift 
compu ted from the 30- band COSMOS photometric cata¬ 
logue dllbert et al.l[200^ . The spectroscopic (photometric) 
sample is complete for group members with M* > 10^®'®M© 
(M* > 10® M©). Again, since most of the contribution to the 
total SHMR originates from M* > 10^®M© galaxies, the 
comparison with our results is fair. We note a slightly lower 
value at high mass, and a good agreement within the error 
bars at the group-scale halo mass. 

The value of the central SHMR peak may also be com¬ 
pared to that of ICoupon et al.l ([2012tl computed from a clus¬ 
tering and galaxy number density analysis of the CFHTLS- 
Wide. In their study the authors have measured the evo¬ 
lution of the SHMR peak as function of redshift and have 
found a lower value compared to ours (1.1 x 10^® M© at red¬ 
shift 2 ~ 0.7). The difference may not be fully explained by 
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Figure 12. Stellar-to-halo mass ratio (SHMR) as function of halo mass compared with observations from the literature. Our best-fit 
result for total (central plus satellites) SHMR is shown as the black shaded area. The black das hed line represents the best-fit centra l 
relationship, whereas the dot-dashed line is for the integrated stellar-mass satellite contribution. For lBehroozi. Wechsler &: Conr^ ll2013l) . 
only the central SHMR was published and we display it here for comparison with ou r central SHMR and as an i llustration of t ypica l 
stellar mass systematics. The length of the grey arrow represents the shift to apply to lLeauthaud et ^ ll2012l ) and iGeorge et aJ] J 2 OIII ) 
to reconcile their results with ours, based on the stellar mass comparison with lllbert et alTil201^ 


cosmic variance, firstly because our field significantly over¬ 
laps with the full CFHTLS and secondly because the dif¬ 
ference is larger than our error bars. In fact, due to their 
selection in the optical {i < 22.5), the SHMR peak above 
2 ; = 0.6 is much less constrained than for our Ks < 22 sam¬ 
ple, and their peak location suffers from higher uncertainties 
than in this study, not properly accounted for in their pub¬ 
lished error bars. 

In Fig. [T3] we compare our results with a number 
of se mi-analytic predictions from t he Millennium simula¬ 
tion llSpringeh Frenk fc Whitell2006l) . In brief, semi-analytic 
models are anchored to the dark matter halo merger trees 
provided by N-body simulations, in which empirical recipes 
of physical processes drive the evolution of galaxies. The 
hne-tuning of those different processes aim at reproducing 
the observed galaxy statistical properties. In each case, to 
derive the total SHMR we compute the sum of the cen¬ 
tral galaxy stellar mass and the integrated stellar masses 
of satellites with M* > 10^^Mq to match our sample mass 
completeness limit. The central SHMR is represented as a 
dashed line and the shaded area represents the total SHMR 
with 15% and 85% percentiles. All qu antities were c ompu ted 
at redshift 2 = 0.8. The model of iBower et al.l ll20 06l) is 
shown in red (top left), the model of iDe Lucia fc_^lmgotJ 
(2007) in orange (top right), and the model of IGuo et al.l 
1 2 OIII) - a mod ified version of the la tter - in green (bottom 
right). In both iDe Lucia fc Blaizotl and IGuo et al.l models, 
the contribution from satellites to the total SHMR is signifi¬ 
cantly below the observations. Despite a different treatment 
of satellite galaxies and the efficiency of stellar feedback in 
the latter model, compared to the former, those changes do 
not show up here. The discrepancy with our results could not 


arise from a limitation caused by the simulation resolution, 
as we imposed a cut of M^, > 10^°M (7) to match our ob¬ 
servations. The model of iBower et al.l better reproduces the 
observed satellite SHMR, however it underestimates the cen¬ 
tral SHMR and features a significant scatter in the M* — Mj, 
relationship. 

We also show the re s ults f rom the analytical model 
proposed by iBirrer et"^ (l2014h in blue (bottom right). 
Their model is a n application of the gas-regulator model 
dmiv et al.ll2013l ). in which the star formation efficiency is 
driven by the amount of available gas in the reservoir. In its 
simplest form, the model describes the inflows and outflows 
of the gas in the reservoir by two adjustable parameters: a 
star-formation efficiency e, and a mass-loading factor A that 
represents the outflows, proportional to the SFR. We show 
their SHMR at 2 = 1 from the model “C”. 


6 DISCUSSION AND CONCLUSIONS 

Using a unique combination of deep optical/NIR data and 
large area, we have combined galaxy clustering, lensing and 
galaxy abundance, to put constraints on the galaxy occupa¬ 
tion function in the range 0.5 < 2 < 1 and to link galaxy 
properties to dark matter halo masses. Our main result is 
an accurate measurement of the central galaxy M* — Mh re¬ 
lationship at 2 ~ 0.8 ranging from halo masses at the peak 
of the SHMR up to the galaxy cluster mass regime. We also 
provide separate measurements of the SHMR for central and 
satellite galaxies. 

We have shown that the statistical errors (computed 
using a jackknife estimator) were smaller than systematic 
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Figure 13. Stellar-to-halo mass ratio (SHMR) as function of halo mass compared with simulations from the literature. We compare 
our total, central and satellite SHMR results with three studies based on semi-analytic models applied to the Millennium simulation 
(top-row and bottom-left panels) and one study (bottom right panel) based on the “gas-regulator” analytical model. For each model we 
also display the corresponding central (dashed line) and satellite (dot-dashed line) SHMR. 


errors in the stellar mass measurements caused by uncertain¬ 
ties in the assumed cosmology, dust modelling, and photo¬ 
metric calibration. Due to the relatively small amount of sta¬ 
tistical uncertainties, the low- to intermediate-mass regime 
of the stellar mass function is most affected by systematic 
errors; a factor of ~ 8 was found between statistical errors 
and total errors, increasing the error bars of parameters con¬ 
trolling the shape of the M* — Mu relationship by approx¬ 
imately the same amount (see Table mi. Conversely, clus¬ 
tering and lensing measurements feature relatively higher 
statistical uncertainties and only a factor of ~ 2 increase in 
error of the HOD parameters describing the satellite popu¬ 
lation is observed compared to statistical errors. By probing 
such a large volume, nearly 0.1 Gpc®, this study brings un¬ 
precedented constraints on the M* — Mh relationship from 
statistical methods in the cluster mass regime at those red- 
shifts. As shown in Fig. [TOl our results make the link be¬ 
tween statistical methods based on HOD applied to deep, 
small-volume surveys, with direct measurements of massive 
clusters from large-scale surveys. 

For central galaxies, we have shown that when properly 
accounting for halo mass definition, choice of the IMF and 
the scatter between M* and Mh, there is a general agree¬ 
ment among results from the literature. We find that stel¬ 
lar mass estimates are the main source of uncertainty, as 
reflect ed by the light-blue shaded area from iBehroozi et al.l 
(120181) in Fig. HOI or the stell ar mass shift measu red between 
iLeauthaud et all ll2012l) and lllbert et al.l (l2010l ). We stress. 


however, that if stellar mass differences may induce a global 
shift (for instance caused by a separate choice for the IMF), 
it may also translate into a mass-dependent shift in the more 
general case (e.g. between two sets of SPS models): hence 
applying a constant shift may not necessarily reconcile two 
measurements. 

In Fig. 111! stellar mass systematics do not seem to ex¬ 
plain all of the observed differences with some results from 
the literature for which the stellar mass was measured in a 
similar way to this study. To measure the impact of some of 
the assumptions made in our model, we have compiled a list 
of potential systematics propagated through the halo mass 
and satellite normalisation best-fit parameters. We quote an 
estimate of 50% error in Mi and 20% error in Baat, respec¬ 
tively. 

For satellite galaxies, the combination of lensing and 
clustering in this work represents a significant improvement 
over studies using only the stellar mass function. In Fig. [12] 
we have shown the measured total SHMR as function of halo 
mass, compared with a number of results from observations 
and simulations in the literature. Starting from group-size 
halos up to the most massive clusters, we find that the to¬ 
tal SHMR is gradually dominated by the contribution from 
satellites. 

Clearly, most SAMs tend to underestimate the total 
amount of stellar mass produced in medium- to high-mass 
satellites < M*/Mq < 10^*^) at 2 ~ 1 compared 

to observations. This would suggest that, in SAMs, the 
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bulk of star formation occurs in low-mass galaxies, but 
is quenched or suppressed at higher mass. Possible expla¬ 
nations for this include either a too strong quenching of 
halos in the mass regime 10^ ° < M*/Mq < 10^^ (e.g. 
the work by iHenriaues et al.l l2012l . who argue that the 
gas could be later reincorporated into the halos), or that 
low-mass sub-halos are too numerous and would “catch” 
the gas in detriment of high-mass sub-halos. It is inter¬ 
esting to link this feature to the overabundance of low- 
mass galaxies found in numerical simu l ations compared 
to observations (see e.g. G_uo_et_^2 20nWeinmann et al.l 


I 2 OI 2 I: [P^ Lucia. Muzzin fc Wein mannl l2014l ~). In this con¬ 
text. ISchive. Chiueh fc Broadhurst (12014 1 recently proposed 
that cold dark matter could behave as a coherent wave and 
have shown using N-body simulations that this would sup¬ 
press a large amount of small-mass halos. 

Finally, we can summarise our findings as follows: 


• the HOD model accurately reproduces the four observ¬ 
ables within the statistical error bars in all mass bins over 
three orders of magnitudes in halo mass and two orders of 
magnitudes in stellar mass; 

• our M* — Mh relationship shows generally good agree¬ 
ment with the literature measurements at 2 ~ 0.8 and we 
have shown that, when modelling differences are properly 
accounted for, we are able to make a fair comparison of a 
number of results derived using independent techniques; 

• the systematic errors affecting our measurements were 
propagated through the whole fitting process. For the pa¬ 
rameters describing the M* — Mh relationship, we find that 
including systematic errors leads to a factor of 8 increase 
in error bars, and for the parameters describing the satel¬ 
lite HOD a factor of 2 increase in error bars, compared to 
statistical error bars; 

• the sum of systematic errors from the halo model and 
our model assumptions may be as high (but likely overesti¬ 
mated) as 50% in halo mass and 20% in the satellite number 
normalisation; 

• the central galaxy stellar-to-halo mass ratio peaks at 
Mh = 1.9 X IO^^Mq, a value slightly l arger than the c luster - 
ingresults from the full CFHTLS from ICoupon et al.l(l2012li . 

• the total (central plus satellites) SHMR is dominated 
by the satellite contribution in the most massive halos, in 
apparent contradiction with SAMs in the Millennium simu¬ 
lation. 


We have demonstrated the power of associating a large 
and deep area with a combination of independent observ¬ 
ables to constrain the galaxy-halo relationship with unprece¬ 
dented accuracy up to 2 = 1. The potential of these data 
will undoubtedly allow us to extend this analysis to galaxies 
split by type in a future work. 

Additionally, studying the evolution in redshift of the 
SHMR above 2 = 1 is one of the greatest challenge in the 
near future. If abundance matching already probes the cen¬ 
tral galaxy-halo relationship up to high redshift, clustering 
and lensing are necessary to put constraints on the satel¬ 
lite HOD and break some of the degeneracies. Large-scale 
clustering measurements require wide-held imaging, whereas 
high-redshift lensing techniques are yet to be improved, but 
on-going projects such as HSC, DES or COSMOS/SPLASH, 
which will increase by orders of magnitude the currently 


available data, represent the ideal data sets to address those 
issues. 
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APPENDIX A: COMPLETENESS OF THE 
SAMPLES 

In this section we use the CFHTLS-Deep/WIRDS combined 
data to test our samples’ mass completeness. The CFHTLS- 
Deep/WIRDS data are over 2 magnitudes deeper in all 
bands compared to our CFHTLS-Wide/WIRCam data and 
with accurate photometric redshift and stellar mass esti¬ 
mates computed in a similar fashion to this study. Fig. ED 
shows the galaxy distribution in WIRDS as function of stel¬ 
lar mass and redshift corresponding to our selection Ks < 22 
for the photometric sample (top) and i < 22.5 for the spec¬ 
troscopic sample (bottom). 

The density fluctuations seen as function of redshift are 
due to cosmic variance (the field of view is smaller than 
1 deg^), but we do not expect any significant impact on our 
completeness assessments. In both panels we represent the 
90% completeness limits as dashed lines, and our samples’ 
selection as red boxes. In the case of the photometric sample, 
a conservative z < 0.7 cut is adopted in the lower mass 
sample to prevent missing red galaxies caused by the optical 
incompleteness at the CFHTLS-Wide depth. Overall, these 
verifications show that all of our samples are complete in 
mass. 



Figure Al. Galaxy distribution as function of stellar mass and 
redshift in WIRDS. Stellar mass 90% completeness limits of Ks < 
22 (top) and i < 22.5 (bottom) selected samples are represented 
as the dashed black line and the sample selection as the thick red 
line. 


APPENDIX B: DETAILS ON THE 
DERIVATION OF THE OBSERVABLES 

Here we provide detailed calculations of the four ob¬ 
servables used in this study and derived from the HOD 
model described in Section 21 For the dark matter halo 
profile and the distribu t ion o f satellites, we assume a 
ikavarro. Frenk fc White! lll997li (NFW) profile with the 
theoretical ma s s-conc entration relation from Eq. (16) of 
iTakada fc JainI ll2003ll with cq = 11 a nd 0 = 0.13, featur- 
ing the redshift dependence (I + z)~^ dBullock et al.lliooifl . 
All dark matter quantities are derived at the mean redshift 
of the galaxy sample, computed from the expectation value 
of the sum of redshift PDFs. All quantities are computed 
in comoving units (“co”). The clustering and galaxy-galaxy 
lensing are then converted into physical units (“phys”) to 
match the measurements. 

B1 Stellar mass function 

The stellar mass function is the integrated HOD over the 
halo mass function: 

</smf(M‘SM^) (B1) 

= ^ (Wot(Mh|M*‘%M**2))^dMh. 


B2 Galaxy clustering 

We describe galaxy clustering using the two-point correla¬ 
tion function, as the sum of the one-halo and two-halo terms: 

^gg(rco) = 1 + ^gg,l(rco) + ^gg,2{rco) ■ (B2) 

The one-halo term, ^gg,i(rco), expresses the relative contri¬ 
bution of galaxy pairs within the halo (Atot(Mh)(Atot(Mh) — 
l))/2 and can be decomposed, assuming Poisson statistics 
for the satellites, into two terms: 

(iVceniVs.t)(Mh) = (Acen(Mh))(Vsat(Mh)); 

(Vsat(A,at-l))(Mh)/2 = (Wat(Mh))V2. (B3) 

The correlation function for central-satellite pairs is 
given by 

1 ^Cs(^C0 5 0 ~ 

[ dMhn(Mh,2)-^^^^^^|^^^Ph(rco|Mh), (B4) 

where we assume that the distribution of central-satellite 
pairs simply follows that of the dark matter halo profile. 
The lower integration limit Mvir(rco) accounts for the fact 
that no halo with a virial mass corresponding to rco would 
contribute to the correlation function. 
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For the satellite contribution ^as, the distribution of 
satellite pairs is the convolution of the dark matter halo pro¬ 
file with itself, computed here in Fourier space. The satellite 
power spectrum is 

PUk)= [ 

■J M,o„ "'gal 

where Uh(fc|Mh) is the Fourier transform of the dark-matter 
halo profile ph{rco\Mh)- The correlation function ^ss is then 
obtained via a Fourier transform. 

The one-halo correlation function is the sum of the two 
contributions, 

^gg,l(rco) ~ 1 -j- ^cs(rco) -b ^ss(rco) . 

The two-halo term is computed from the galaxy power 
spectrum: 


T2(fc,rco) = T’m(fc) X 


(B7) 


(^co) 


(iVtot) 


dMhn(Mh)-7——r&h(Mh,rco)|uh(fe|Mh)| 


L-'Mio 


‘gal 


(r. 


where 


r*flim(’'co) 

’^gal(^co) = / n{Mh){Ntot} dMh- 

Jmi„ 


(B8) 


The upper integration limit Mumfrco) acco unts for halo ex¬ 
clusion as detailed in ICoupon et al.l (|2ni2ll , and references 
therein. 

Finally, the two-halo term ^gg ,2 of the galaxy auto¬ 
correlation function is the Fourier transform of Eq. (lB7l) 
renormalised to the total number of galaxy pairs: 


1 + Cgg.2(l'co) 


n'gaiirco) 

^gal 


[1 + Cgg,2(’'co)] . 


(B9) 


The projected clustering w{8) is derived from the pro¬ 
jection of ^gg onto the estimated redshift distribution from 
the sum of PDFs, assuming t he Limber approximation (see 
details in lCoiiPon et al.ll2012l L 

The real-space clustering Wp(rp,co) is derived from the 
projection of the 3D correlation function along the line of 
sight: 


UJp(^p,co) — 2 



rcodrcoCgg(rco) (rco - rp,co) , (BIO) 


converted into physical units as: 


'^^p,phys — rCp^co/(1 ^) • 


(Bll) 


B3 Galaxy-galaxy lensing 

The galaxy-galaxy lensing estimator measures the excess 
surface density of the projected dark matter halo profile: 

^Eco(r'p,co) ~ Eco(^ r*p,co) Eco(r*p,co) , (7112) 

where Eco(< rp,co) is the projected mean surface density 
within the comoving radius rp,co and Eco(rp,co) the mean 
surface density at the radius Cp^co. 

To compute the analytical projected dark matter den¬ 
sity E, we write 



— P 1 “b ^gm ^ 


dTTco, (B13) 


where rp,co is the transverse comoving distance, ttco the line- 
of-sight comoving distance, p the mean density of the Uni¬ 
verse, so that AEco(rp,co) is related to the galaxy-dark mat¬ 
ter cross-correlation function ^gm through 


AEco(rpiCo) — Eco(^ ^p,co) Eco(r’p^co) — 

^ pPp,co p'^Tn&x. 


(B14) 


^ pPp,co P'^Tn&x , _^ 

2 / / ^p,co^gm f ■\/^p,co H“ '^c^j dTTcodrp (. 

P'^m&x / _ 

“2 J ^gm ^ 's/ ^p,co TTco J dTTco 


The integration along the line of sight is performed up to 
the scale TTmax = 80 Mpc. 

The excess surface density in physical units writes: 


^Sphys — ASco X (1 -}- Zl^) , 


(B15) 


where zl is the redshift of the lens galaxy. 

As for ^gg, ^gm can be written as the sum of the one- 
and two-halo terms: 


Cgm(r) — 1 -b 5gm,i(r) -b ^gm, 2 (r). (B16) 

Cgm,i(r) is itself decomposed into a contribution from the 
cross correlation of the central galaxy-dark matter and from 
that of the satellite-dark matter, both assuming a NFW pro¬ 
file. We write the former as 


1 + Cgm,cen(r, z) — 


f'^high 




dMhn(Mh,z)A^pj^(p|Mh)^ (B17) 
^gal P 


and the latter (^gm,sat from the Fourier transform of its power 
spectrum 


Fgm,ss (^) — 

f dMhn(Mh) ^\uh{k\Mi, f . (B18) 


Finally we compute the two-halo term ^gm, 2 (r) from the 
Fourier transform of the galaxy-dark matter cross correla¬ 
tion power spectrum: 


T’gm,2(fc,r) = Pni(fc) (B19) 

X / dMhn(Mh)A2ii^bj^(Mh,r)|uh(fc|Mh)|, 

"gal(0 

with a similar treatment of halo exclusion to that of the 
galaxy power spectrum. 


APPENDIX C: SYSTEMATICS CHECKS ON 
LENSING AND CLUSTERING 

We have performed systematics checks for the lensing 
and clustering measurements. In Fig. ET] we detail the 
galaxy-galaxy lensing measurement for the sample 10.40 < 
\og{Mi,/M q) < 10.65 as an example. The top panel shows 
the data (dots with error bars) and best-fit model (thick 
line) with the different components of the model the cen¬ 
tral galaxy term, the satellite term, and the 2-halo term. 
The lower panels show a number of systematics checks. The 
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Figure C2. Galaxy-galaxy lensing measurements separating the 
background sample into 0.8 < 2 p < 1.2 sources (purple dots) and 
2 p > 1.2 sources (green triangles), keeping the same lens galaxy 
foreground sample (low-redshift galaxies with spectroscopic red- 
shifts). 



Figure Cl. Galaxy-galaxy lensing measurements and system- 
atics checks for the sample 10.40 < Mq) < 10.65. In 

the top panel we show the data (dots with error bars) and the 
model (thick line) split into the stellar term in dotted line, the 
central term in dashed line, the satellite term in dot-dashed line 
and the 2-halo term in triple dot-dashed line. The lower panels 
show the systematic tests (rotated-shape signal and random lens 
positions), calibration factor (multiplicative bias correction and 
boost factor), and the lower left corner the correlation coefficients 
of the correlation matrix from the jackknife estimate. 


“ex” panel shows the signal measured after rotating the el- 
lipticities by 45° and the “ran. lenses” panel shows the signal 
measured by randomising the lenses positions, both consis¬ 
tent with zero. The “H-m” panel shows the multiplicative 
bias correction applied to the galaxy-galaxy lensing mea¬ 
surement, estimated after replacing the ellipticities by the 
multiplicative calibration factor 1 -|- m. The “boost factor” 
was estimated from randomising the background source po¬ 
sitions and measuring the ratio of the number of real sources 
over random objects as a function of distance from the 
lenses, and applied to the galaxy-galaxy lensing measure¬ 
ment. The covariance matrix from the jackknife estimate is 
shown in the left bottom corner of the figure. The relatively 
small off-diagonal values show the low correlation between 
data points. We repeated identical tests for all mass bins. In 
all cases, systematics are found to be consistent with zero. 

In Fig. IC2l we test the impact of including high redshift 
sources beyond 2 > 1.2. To do so, we select an arbitrary sam¬ 
ple of low redshift lens galaxies with a spectroscopic redshift 


and we measured the galaxy-galaxy lensing signal using all 
sources with 0.8 < 2 p < 1.2 (purple dots in the figure) and 
all sources with 2 p > 1.2 (green triangles in the hgure). We 
see no significant difference between the two signals, mean¬ 
ing that the photometric redshifts and shape measurements 
in our catalogue are robust enough beyond 2 p > 1.2. 

In Fig. lC3l we show the projected clustering in the mass 
bin 10.60 < \og{Mt/M q) < 10.80. The top panel shows the 
data points with error bars and the best-fit model, with the 
different components of the model: the 1-halo term split into 
the central-satellite and satellite-satellite terms, and the 2- 
halo term. In the middle panel we show the corresponding 
HOD, as a dashed line for the central contribution and as a 
dot-dashed line for the satellites’ contribution. 


APPENDIX D: 2D CONTOURS 

We show in Fig. IDll the likelihood distributions of the best- 
fit HOD parameters. Here the results are shown for the 
MCMC run done with total (statistical plus systematic) er¬ 
rors. 
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Figure Dl. ID (diagonal) and 2D likelihood distributions of best-fit HOD parameters in the case of total errors. The 2D contours 
represent the 68.3%, 95.5% and 99.7% confidence limits. We used flat priors within the ranges shown on the figure for all parameters. 
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Figure C3. wiO) measurements and the corresponding HOD 
function for the sample 10.60 < log(M*/M 0 ) < 10.80. In the 
top panel, we show the data points with error bars and the best- 
fit model: the dotted line represents the central-satellite cross¬ 
correlation, the dashed line the satellite-satellite auto-correlation, 
and the dot-dashed line the central-central auto-correlation (or 2 - 
halo term). The middle panel displays the corresponding HOD, 
the dashed line shows the central galaxy HOD and the dot-dashed 
line the satellites’ HOD. The lower right panel shows the corre¬ 
sponding redshift distribution constructed from the sum of indi¬ 
vidual PDFs. The lower left panel shows the correlation coeffi¬ 
cients of the covariance matrix from the jackknife estimate. 
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